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1.  Introduction. 

In  the  field  of  Fracture  Mechanics  not  very  much  theoretical  work 
has  been  done  in  order  to  assess  analytically  the  three-dimensional 
stress  character  which  prevails  at  the  base  of  a  stationary  crack.  As 
a  result,  most  of  our  current  design  criteria  are  based  on  already 
existing  two-dimensional  solutions  and  therefore  are  in  general 
inadequate.  For  example,  the  common  experimental  observation  of  a 
change  from  ductile  failure  at  the  edge  to  brittle  fracture  at  the  center 
of  a  broken  sheet  material  has  so  far  defied  analysis.  Yet  an  orderly 
theoretical  attack  on  the  problem  can  provide  important  guidance  to 
this  and  other  phases  of  fracture  research. 

Hie  mathematical  difficulties,  however,  posed  by  three-dimensional 
fracture  problems  are  substantially  greater  than  those  associated 
with  plane  stress  or  plane  strain.  Be  that  as  it  may,  the  author  would 
like  to  investigate  the  subject  further  at  least  within  the  theory  of 
linear  elasticity.  While  he  recognizes  the  fact  that  this  theory  can¬ 
not  include  the  nonelastic  behavior  of  the  material  at  the  crack  tip 
per  se,  it  can  evince  many  characteristics  of  the  actual  behavior  of 
a  cracked  plate,  including  those  due  to  thickness.  Thus  the  theory  of 
elasticity  is  a  logical  fountainhead  for  detailed  theoretical  study. 


2.  Historical  Develi 


opment. 

There  exist  in  the  literature  very  few  analytical  papers  that  deal 
specifically  with  the  three-dimensional  stress  character  at  the  base  of 
a  stationary  crack.  Moreover,  in  their  present  form  these  papers  are 
not  only  incomplete  but  also  contradictory.  As  a  result,  much  con¬ 
troversy  and  many  doubts  have  been  raised.  It  is  appropriate,  therefore, 
to  discuss  these  papers  and  their  respective  results  in  chronological 
order. 

In  1972,  Benthem,  using  the  method  of  separation  of  variables*, 
was  able  to  solve  for  the  stress  distribution  in  the  neighborhood  of 
the  comer  point**  of  a  quarter  plane  crack.  His  results  [2]  show  that 
the  stresses  there  behave  like  p"a  ,  where  0.500  <  a  <  0.709.  In  order  to 
obtain  the  order  of  the  singularity,  Benthem  had  to  trancate  an  infinite 
system  which,  in  turn,  he  solved  for  the  eigenvalues  numerically.  This 
approach,  however,  raises  three  important  questions:  One,  is  the 
solution  really  separable,  particularly  in  0  and  <p  ?  Two,  is  the 
solution  thus  obtained  complete?  Three,  should  the  numerical  deter¬ 
mination  of  the  singularity  from  a  truncated  system  be  trusted? 
Unfortunately,  Benthem  has  provided  no  answers  to  any  of  the  above 
important  and  difficult  questions. 

A  few  years  later,  Folias,  using  a  method  developed  by  Lur'e  [3] 
apd.  the  application  of  Fourier  Integral  Transforms,  was  able  to  solve 
[4]  Navier's  equations  for  a  more  complicated  problem,  that  of  the 

*This  method  was  fully  articulated  by  M.L.  Williams  [1]  for  classical 
planar  elasticity* in  order  to  establish  the  singular  behavior  at 
re-entrant  comers. 

**That  is  the  point  where  the  crack  front  meets  the  free  surface  of 
the  half  space. 
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3-D  Griffith  crack  (see  Figure  1).  The  integrals  were  subsequently 
expanded  asymptotically  and  the  stress  field,  valid  in  the  very  inner 
layers*  of  the  plate,  was  recovered.  From  the  results,  one  concludes 
that  in  the  very  inner  layers  of  the  plate: 

(1)  the  stresses  possess  the  usual  singularity, 

(2)  the  stresses  posses  the  usual  angular  distribution, 

(3)  the  stress  intensity  factor  Kj  is  a  function  of  z  , 

(4)  exact  plane  strain  conditions  exist  only  on  the  plane  z  ■  0  , 

(5)  a  pseudo  plane  strain  state  exists  and  the  equation 

az  *  v(>x  +  V 

is  satisfied, 

(6)  as  the  plate  thickness  2h  -*■  <*>  ,  the  plane  solution  is  recovered, 

(7)  as  Poisson's' ratio  v  -*•  0  ,  the  plane  stress  solution  is  recovered. 
Furthermore,  he  was  able  to  show  that  at  the  corner  the  stresses  are 
proportional  to 


p-(H  +  2v) 


(6,0). 


In  order  to  recover  the  value  of  the  singularity,  Folias  solved  analyti 
cally  a  difference-differential  equation.  Unfortunately,  because  of 
the  enormous  difficulties  which  the  integral  representations  presented 
at  the  comer,  he  was  unable  at  the  time  to  recover  the  functions 
fij(6>40  explicitly. 


*The  reader  should  note  that  the  asymptotic  expansions  are  only  valid 
for  (z/h)  «  1  and  for  c/h  «  1.  This  is  because  h  was  assumed 
to  be  very  large  so  that  a  perturbation  about  the  well-known  plane- 
strain  solution  could  be  made. 
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It  should  be  emphasized  that  Folias's  main  result  at  the  comer 
should  be  interpreted  as  "the  singularity  at  the  comer  can  at  most  be 
of  the  order  (%  +  2v)"  .  This  is  because  the  functions  f  —  (0,<tO 
could  very  well  be  of  the  type  that  do  vanish*  in  th^  neighborhood  of 
the  comer  point.  Thus  Folias's  result  may  or  may  not  be  in  contra¬ 
diction  with  Benthem's. 

Researchers  in  the  field  of  Fracture  Mechanics,  however,  were 
unwilling  to  accept  the  possibility  of  an  infinite  displacement  field 
on  the  basis  of  physical  intuition.  Consequently,  the  results  were 
considered  highly  controversial  and  the  following  two  legitimate  questions 
were  raised**:  Is  the  solution  really  complete?  Two,  do  the  series 
representations  converge?  Unfortunately,  Folias  provided  no  answers 
to  any  of  the  above  questions. 

In  1976,  Kawai  [7],  using  the  method  of  separation  of  variables 
was  able  to  obtain  an  alternate  solution  to  Benthem's  problem.  Although 
the  method  of  approach  is  essentially  the  same  as  that  of  Benthem's, 
his  results  are  definitely  contradictory***.  His  results  show  that  at 
the  comer  the  stresses  behave  like  p’a  ,  where  \  <  a  <  1  .  In 
determining  the  singularity,  Kawai  used  the  collocation  method  in  order 
to  satisfy  the  three  boundary  conditon  on  the  free  surface.  Thus,  as 
in  Benthem's  case,  the  same  questions  apply  to  this  work  also. 


*The  reader  should  note  that  this  result  was  actually  obtained  by 
'marching  out'  the  solution  from  the  inner  to  the  outer  layers,  and  as 
a  result  such  a  hypothesis  may  not  be  totally  unreasonable.  See  also 
comments  on  p.  5. 

**See  Discussion  of  paper  by  Benthem  and  Koiter  [5]  and  author's  Closure 

[6]. 

***Ma thematically,  Kawai 's  method  of  construction  of  the  solution  is 
more  systematic  than  that  of  Benthem's. 


A  few  months  later,  Benthem discovered  that  his  previously  reported 
solution  was  incomplete  and  that  his  new  results  [8]  now  read 

°ij  ~  with  0  <  a  <  55  ‘ 

Here  again,  the  same  questions  raised  during  his  previous  work  apply 
too. 

Finally,  in  1977  Kawai  [9]  reported  an  error  in  his  previous 
analysis  and  although  the  correction  affected  slightly  the  value  of 
a  the  trend  essentially  remained  the  same. 

In  the  meantime,  Folias  also  discovered  that  his  solution  of  the 
difference-differential  equation  was  not  quite  complete  either*.  The 
correction,  however,  does  not  directly  alter  the  basic  result  at  the 
comer. 

It  is  interesting  to  note  that  Kawai  does  recover  the  same  sin¬ 
gularity  that  Folias  reported.  The  singularity  O  *r-A>),  however,  disappears 
as  he  considers  more  and  more  terms  in  his  collocation  scheme  but  at 
the  same  time  he  experiences  convergence  problems.  Ihis  observation 
strengthens,  perhaps,  the  hypothesis  that  Folias’ s  f —  (9,4>)  functions 
do  indeed  vanish  in  the  neighborhood  of  the  comer  point  and  that  most 
likely  are  needed  in  the  very  inner  layers  of  the  plate.  The  later  has 
also  been  observed  by  Newman  [10]  for  (c/h)  ratios  less  than  one, 
which  is  comparable  with  the  asymptotic  expansion  used  by  Folias. 

*This  is  not  to  be  confused  with  the  question  of  completeness  of 
the  solution  to  Navier’s  equations,  i.e.  eqs.  (52) -(54)  ref.  [4], 

The  corrected  result  to  eq.  (85)  of  reference  [4]  is  given  in  Appendix  j. 
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Be  that  as  it  may,  the  presence  of  a  third  solution  obscurred  the 
issue  even  further  and  essentially  raised  more  questions  than  gave 
answers.  So  the  controversy  still  remains. 


3.  Purpose  of  Present  Work. 

In  view  of  the  preceding,  it  is  evident  that  mathematical  rigour 
becomes  essential  if  one  is  to  avoid  any  possible  pitfalls.  As  a 
result,  the  author  decided  to  seek  the  answers  to  the  following  two 
important  questions  first: 

(i)  Is  the  solution  of  this  notoriously  difficult  problem 
unique?  And  if  so,  under  what  conditions? 

(ii)  Is  the  solution  to  Navier’s  equations  as  given  by  the 
author  in  reference  [4],  i.e.  eqs.  (52)- (54),  general 
enough  to  represent  the  solution  of  this  practical 
problem? 

The  answers  to  both  of  the  above  questions  were  given  by  Prof.  Calvin 
Wilcox. 

* 

First  of  all,  he  was  successful  in  proving  [11]  that  a  displace¬ 
ment  field  that  satisfies  the  condition  of  local  finite  energy  is 
unique.  This  of  course  is  quite  a  departure  from  our  traditional  ! 
2-D  fracture  mechanics  thinking,  for  the  displacements  now  can  be  , 
allowed  to  be  singular.  Consequently,  one  may  not  apriori  assume  them 
to  be  finite  as  it  is  customerily  done.  In  general,  such  an  assumption 
makes  the  class  of  solutions  too  restrictive  and,  as  a  result,  one  may 
not  find  a  solution  to  the  problem.  On  the  other  hand,  the  solution 
could  very  well  give  finite  displacements  everywhere!  Be  that  as  it 
may,  physical  intuition  should  be  used  with  extreme  caution. 

Second,  he  was  able  to  show  [12]  that  the  Fourier  integral  expres- 

JL  ^ 

sions*  representing  the  general  solution  to  Navier’s  equations  are 

♦  S*-«-  oJLSo  3EC  . 

*  #"*See  equations  (52) -(54)  of  reference  [4]. 

o  jH.  • 


**  V 
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complete  and,  furthermore,  the  'symbolic  method'  used  is  justifiable. 
In  order  to  prove  this,  he  used  a  double  Fourier  integral  transform  in 
x  and  y  and  subsequently  a  contour  integration  to  recover  precisely 
the  same  expressions  as  those  reported  by  Folias  in  reference  [4], 

Finally,  it  remains  to  determine  explicitly  the  stress  field 
ahead  of  the  crack  tip  and  throughout  the  thickness  of  the  plate. 

In  reference  [4],  the  author,  by  the  use  of  analytic  continuation, 
attempted  to  'march  out'  the  solution  from  the  inner  to  the  outer 
layers  of  the  plate.  Although  in  principle  this  seems  feasible,  in 
practice  it  is  very  difficult  and  most  of  all  tedious.  Moreover, 
questions  of  convergence  will  inevitably  be  raised.  As  a  result, 
in  this  paper  we  will  use  an  alternate  and  more  elegant  approach  in 
order  to  complete  the  problem. 

By  finding  the  biorthogonal  relation  for  the  eigenvectors,  we 
will  set  up  a  double  integral  equation  for  the  unknown  function  v  , 
which,  physically,  represents  the  projection  of  the  displacement  v 

onto  the  xz-plane.  The  advantages  of  this  new  approached  over  that 
of  reference  [4]  are: 

(i)  we  are  now  seeking  the  solution  to  one  equation  only, 

(ii)  the  unknown  function  is  real  and  furthermore  has  physical 
meaning, 

(iii)  the  kernel  of  the  integral  equation  is  independent  of  the 
shape  of  the  crack*. 


*In  this  analysis  we  restrict  ourselves  to  planar  and  symmetric  cracks 
subjected  to  mode  I  loadings. 


4.  Formulation  of  the  Problem. 

Consider  the  equilibrium  of  a  homogeneous,  isotropic,  elastic 
plate  which  occupies  the  space  |x|  <  ®  ,  |y|  <  ®  ,  |z|  <  h  and 
contains  a  plane  crack  in  the  x-z-  plane  (see  Figure  1).  The  crack 
faces,  defined  by  |x|  <  c  ,  y  *  (T  ,  z  <  h  ,  and  the  plate  faces 
| z |  *  h  are  free  of  stress  and  constraint.  Loading  is  applied  on  the 
periphery  of  the  plate  |x|  ,  |y|  -*■  «  and  is  given  by 

a  *  r  “  x  *  0  .  a  *  . 

x  xy  yz  *  y  °0 

In  the  absence  of  body  forces,  the  coupled  differential  equations 
governing  the  displacement  functions  u,  v,  and  w  are 

•  w  -  ■&>'  *  72tu'v-”)  ■ 0 

2 

where  V  is  the  Laplacian  operator,  m  =  1/v  ,  v  is  Poisson's  ratio. 


e  =  2ii  +  23L  ♦ 
e"  3x  3y  3z 


(4) 


and  the  stress-displacement  relations  are  given  by  Hooke's  law  as: 

(5) -(10) 


«x-2* 


•  ' 

3u  .  e 

1 

1 

5 

07 

c 

<  . 

5x  +  mT 

>  - 

f  »  Txy  ^3y  *  lx] 

with  G  being  the  shear  modulus. 
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-4.  <1 


As  to  boundary  conditions,  one  must  require  that  at: 


1*1 

<  c  ,  y  -  0*  ,  1 z J  <  h: 

T  *  T  ■0*0 
xy  yz  y 

(11) 

1*1 

-  h  : 

T  ■  T  ■  0  *  0 

xz  yz  z 

(12) 

lyl 

■+•  «  and  all  x  : 

T  *  T  *0.0  ■  o"- 
xy  yz  *  y  0 

(13) 

1*1 

eo  : 

0  ■  T  *  T  *0. 

x  xy  zx 

(14) 

It  is  found  convenient  to  seek  the  solution  to  the  crack  plate 
problem  in  the  form 


u  *  u 


CP)  ♦  U(C) 


etc., 


(15) 


where  the  first  component  represents  the  usual  "undisturbed"  or 
"particular"  solution  of  a  plate  without  the  presence  of  a  crack.  Such 
a  particular  solution  can  easily  be  constructed  and  for  the  particular 
problem  at  hand  is 


“tP)  -  -  OT  <"-2>2  *  • 


(pi  2  °nCra’2) 

VCp)  -  -  [1  -  (m-lH-^-y 


w(p)  -  -  (m-2)‘^ 


2*0. 


(16) 


A  -  (m-1)3  -  3(m-l)  ♦  2  . 


where 
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5.  Method  of  Solution. 

The  complementary  solution  to  Navier's  equations,  subject  to  the 
corresponding  boundary  conditions  (12)  and  (13),  is  given  by  reference 
[4]  as: 

(i)  coeglementary  displacements* 

U<C)  ‘  i>!  *  |y|  91  *  PT  e'S,yl 


1  r„  ~ 


-/sV|y| 


517  v£i‘» 


cos ($vh) [ (m- 2+m  cos4 (Bvh) )  cos (Bvz) 


(17) 


-  mB  z  sin (8  z)]  ♦  l  S  e 


n-1 


cos (anz)}sin(xs)ds 


r<c)  ■  T  -  P1  -  )ylQi  -  FT  ”2V  '‘s|yl 


•  -/s  *8^ 

*  ^-rv  * - - - cos(Bvh)  [(ra-2  ♦  m  cos2(8vh))cos(8vz)  - 


(18) 


5PT  v^l  v 


-i4Z+a2|y| 


mB  z  sin(B  z)]  -  J  S  n — Te  cos  (a  z)}  cos(xs)ds 

n-1  /s  ♦<* 


■w  ■  [«&  *  V  '•5|yl*^jlVv 


cos  (B  h) 


(19) 

•  [(2m  -  2  -  m  cos 2 (8  h))  sin(B  z)  -  mB  z  costtz)]}  cos(xs)ds 


•These  complementary  displacements  represent  a  'general  enough',  or 
'complete',  solution  for  the  satisfaction  of  the  remaining  boundary 
conditions.  For  a  discussion  of  this,  see  reference  [12]. 
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with  corresponding  stress: 
(ii)  complementary  stresses 


<L(c) 

-%rm 


cos(6vh)[sin2(6vh)cos(6vz)+evz  sin(Byz)]  • 
•  cos(xs)ds 


(20) 


Txz(c) 

—TT~ 


(21) 


t  (O 

G 


j:* 


™  ii 


Br: 

V  V 


cos (8yh) [cos 


(Bvh)sin(8yz) 


(22) 


♦  8yz  cos (8yz) ] 


m 

1 

n«l 


sin(anz))  cos(xs)ds. 


* 


•  frtsPi  *  lyi*Qi  *  PT  *2l\  -  Sir  <y  «'s|y| 


♦  =2- j,  s2r 


m-2  v*l  pv  v 


s/s  >6 


cos (8  h)  cos (62) 


m-2  v*l 


cos(8  h)[(m-2+m  cos*(B  h)  cos  (8  2) 


-m8v2  sin(6v2)]  . 


"  c 

*  Ji s- se 


cos(an2)>  cos(xs)ds 


JTO”“  f  {('^Q1  '  SP1  '  ,y|sQl  '  nFT  e’S|xl 


+  kJi  Vv ?  cos(8vh)  C0S(svz) 


F7.1  rv 


*5?  -^sTi/i 


_v  e  "V 


cos(8vh)[(m-2+m  cos  (6vh))  cos(8vz) 


-  m6vz  sin(6vz)] 


-/sV|y| 


"  I  S  s  e  n  cos(o z)}  cos(xs)ds 


^  ■ 7  j!  «2<!siNi  *  2sPi  *  pr  z**\  *  2M*V  e's|yl 


2  r  -^sjlyi  2 

-  5T7  I  rv  e  cos(8vh)  [ (m-2+m  cos^(8vh))  cos(&vz) 


v-1 


-  mevz  sin(8vz)] 


C25) 


•  2s ‘‘♦cl'  -/sZ^|y| 


2.  2 


*  Ji 


e  n  cos(anz)}  sin(xs)ds  , 


where  the  ±  signs  refer  to  y  >  0  and  y  <  0  respectively  and  the 

constants  P,  ,  Q.  ,  r  and  S  are  to  be  determined  from  the  remain- 
11  v  n 

ing  boundary  conditions.  Moreover,  an  “  X  (n*l,2,3...),  and  8v 
are  the  roots  of  the  equation 

sin(2Bvh)  -  -(2Svh)  '.  (26) 

This  equation  has  an  infinite  number  of  complex  roots  which  appear  in 
groups  of  four,  one  in  each  quadrant  of  the  complex  plane  and  only  two 
of  each  group  of  four  roots  are  relevant  to  the  present  work.  These 
are  chosen  to  be  the  complex  conjugate  pairs  with  positive  real  parts. 
The  only  real  root  6y  ■  0  must  be  ignored*. 


*The  first  few  roots  are  tabulated  in  Appendix  II. 


'  .♦ 


-  &  Wi 


15 

By  direct  substitution,  it  can  easily  be  ascertained  that  the 
above  complementary  displacements  satisfy  Navier's  equations  and 

furthermore  the  corresponding  stresses  ,  t  ^  ,  t  do  vanish 

z  Za  yz 

at  the  plate  faces  z  -  ±h  . 

Finally,  if  we  consider  the  following  two  combinations  to  vanish  * 


ao  p 

I  ~-cos(evh)[  sin2(Bvh)  cos(3vz)  ♦  sin(Bvz)] 

V*1 


m  r 


(27) 


-  0 


and 


=T y  l  r  cos (B  h)  [  (m-2*m  cos2 (8  h))  cos(B  z)  -  mB  z  sin(B  z)] 
v-1 


(28) 


»  2s2+a2 


l  -jnr$  sn  cosCV)  ‘  ra  ‘  2sPi  ' 2  Sr  qi  "  ° 

n-l  »s  +ct 
n 


m-1 


for  all  |z|  <  h  ,  then  two  of  the  remaining  stress  boundary  conditions 
are  satisfied  automatically,  i.e. 


■  t  -  0  for  all  x  ,  |z|  <  h  and  y  -  0  . 

xy  yz 


We  will  suppress  for  the  time  being  the  satisfaction  of  the  last 
boundary  condition  and  will  focus  our  attention  to  the  continuity 
conditions . 

As  it  can  easily  be  seen,  all  continuity  conditons  are  satisfied 
if  one  considers  the  following  two  combinations  to  vanish 

•^Notice  that  the  derivative  of  eq.  (27)  with  respect  to  z  leads  to  the 
integrand  of  eq.  (22). 
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3**^.  -  p  1  sz20 

L  1  ¥TT  *i  5+T  sz  Si 


w  r. 


(29) 


+  srr  l  cos(B  h)[(2m-2+m  sin2(B  h))cos(B, ,z) 


v*l 


.  »  r 

+  mByz  sin(Bvz)]  +  ~  }  cos(xs)ds  *  0  ; 


I  x  |  >  c  ,V  I  z  I  <  h  . 


and 


f pH  *  Jr  j.  7  6v  i»  *  ^2(SVW)  “»(V) 

A  V*1  S 


(30) 


-  Byz  sin(ByZ)]}  sin(xs)ds  *  0  ;  |x|  >  c  ,v  |z|  <  h  . 


which  by  Fourier  inversion  lead  to: 


1ST  -T  "  pi  "  S7T  «2Qi  +  5TT  l  —  cos(S  h)  • 
m+1  s  1  m+1  u.  m-Z  s  v  * 


•  [(2m-2+m  sin  (Byh))  cos(Svz)  +  mBvz  sin(Byz)]  (31) 


4m 

+  -—■» 
ra-2 


l  "r*  r  v(^»°»z)  cos (sC)d^ 

V»1  J  "  JQ 


and  * 


*The  reader  should  note  that  eqs  (31)  and  (32)  automatically  satisfy 
eq.  (28). 


w> 


(32) 


r 

'  5+TQl  +  il?  \  “7  ev  cos(evh)[Cl  +  cos2(evh))cos(evz) 

V*1  s 

-  Sv2  sin(6u2)]  ■  7  I  }0  (§7  -  IfVo  s“(s«)d« 

'  4  f  f'  2#y-0  s“ts«« 

“  ~  1  v(5»°>2)  cos(s5)dC  . 

*  Jo 

Adopting  next  the  following  definitions 

Z^1)(z)=(3vh)2[8vh  sin(Svh)  cos(@vz)  -  Bvz  cos(Bv h)  sin(8vz)]  (33) 

Z^3J(z)  =-(8vh)2[8uh  sin(8yh)  cos(Byz)  -  Bvz  cos(Bvh)  sin(Bvz)]  (34) 
-  2(Bvh)2cos(Bvh)  cos(Bvz) 

f (1)  v  5  T  §2<rr>  |„  '  E s2v}  cos  (s»d5  (35) 

2  eC. 

f(3)(z)  i  ±  —  (^  |q  v  cosCs^  -  2&  Qih2  ’  C36) 

equation  (32)  and  the  second  derivative  of  equation  (31)  with  respect 
to  z  become 

2 

•jjrl  7  zfJ(z)  -  f°  v«,0,  z)  cos(sC)dC  (37) 

m  v»l  s  v  ^0 
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and* 

'  ^  X7{z‘1>w 

' 1  ^  (I7 '  sH-« cosCs5)d5  (38) 

*  2#)  »lh2 


respectively,  which  upon  simplifying  one  has 


«  r 

I  (4) 

vl  s 


fC3) (Z) 

.  f(1)(z) . 

(39) 


Next,  following  reference  [13],  we  can  construct  the  biorthognal 
relations 


w^fz)  =  -  g*z  cos  (g*h)  sin(g*z)  +  g*h  sin(g*h)  cos(g*z)  (40) 


and 

W^(z)  =  -  S*z  cos (g*h)  sin(g*z)  +  [g*h  sin(g*h)  (41) 

-  2  cos (g*h) ]  cos  (g*z)  , 


where  g*  stands  for  the  complex  conjugate  of  the  g^  roots.  The 
orthogonality  condition  now  reads 

♦Notice  that  the  continuity  conditions  are  to  be  satisfied  in  the 
interior  of  the  plate  only. 


i 


or 


(jr)  eKv  =  K  fh{WvC4)cn)  £C1)(n)  +  [2  w*C4)(n)  ‘ 


(43) 


f^(n)}dn 


where  for  simplicity  we  have  defined 


/v  2  E 


-h 


(w*(4)(n)  z^cn)  *  [2w*w(n)  -  w*(2)(n)] 


(44) 


.  z<3) 


(n)  )dn  . 


Finally,  in  view  of  equations  (33) -(36),  (40) -(41)  and  (43) -(44), 
one  finds  after  some  simple  calculations  that 


&•-*  oyo2  =os4(6vw 


(45) 


and 


I 


jkv  Hf)  “ 4  he  i:  |.hvC5,0’n)  cos(s° 


{(hs)  ievh  sin(Svh)  cos(Bvn)  -  Byn  cos(Bvh)  sin(Bvn) 


+  2  (-2p)  cos(Bvh)  cos(Bvn)l 


+  0*v)  [0vh  sin(Svh)  cos(Byn)  -  Byn  cos(Bvh)  sin(Bvn) 


+  2  cos(Bvh)  cos(Bvn)]>  dn  d£  . 


Similarly,  from  equations  (27)  and  (32),  we  find  that 


-Ja—(anh)2  f  fh  v(S,0,n)cos(s5)  cos(ann)  dn  d?  (47a) 

/s4+a  '0  '-h 


Ql  ■  *  -’gf-  |  |  v(C,0,n)  cos(sC)  dn  d£  . 


Returning  now  to  the  last  boundary  condition,  we  require  that* 


*Where  we  have  made  use  of  eq.  (28). 
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l0l  'Q1  7'^^'CMCB''WU1  *  “’V”  “S(6V° 


-  Bvz  sin(evz)] 


,  »  r  3  - 

+  — r  y  —7  [  j  n 


-  m  Bvz  sin(Bvz)] 


-  s  ]  cos(Bvh) [(m-2  +  m  cos  (Bvh))  cosfB^z) 


'(48) 


2 

% 

T3 


which,  upon  lasing  the  relations  (46) - (47)  and  interchanging  the  order 
of  integration,  can  also  be  written  in  the  form  of  a  double  integral 
equation  i.e., 


^  1 1  |  f  {±  v^'°’n}  a|5i-H1[|x.-5|;n,z]  dnd£ 
V=1  cracK 
faces 


faces 


(±  v(C,0,n)}  H2[ |x-c|  ;n,z]  dnd? 


-4|  f  U  V(5,0,n))(^)  ^  [^IdndC 

crack 

faces 


(49) 


) 


* 


0  . 
2U’ 


| z |  <  h  ,  |x|  <  c 


•where 


m  B  h 

H,  [  |x-5|  ;n,z]  -  A-  ~y~  [cos (3  h)  sin($vz) 

e  v 

+  8vh  sin(8vh)  sin(6vz)  +  3yz  cos(Bvh)  cos(Bvz)]  • 

•  [ev|x-e|  ]  [evh  sin(evh)  cos(Bvn)  -  evn  cos(Bvh)  sin(Bvn) 

♦  2  (Sjji)  cos(8vh)  cos(3vn)) 


r  5  Ko[0jx'|  ]  dx*  •  [8vh  sin(8vh)  cos(Bvn) 

Jo 


-  cos(B^h)  sin(8vn)  +  2  cosfB^h)  cos(Bvn)]> 


+  5PT  fTHT  f2! T  cos(8vh)  sin(8v2)  *  8vh  sin(3vh) 

ve  v 

&y 

sin(Bvz)  +  Bvz  cos(Bvh)  cos(Bvz)]  •  {[  -  jjqr  Kg  [Bjx-Sl  ] 


'  evh2|jpffKilBJx‘5lI  ■  Ix-^-TF  (6Jx'5l) 


2B,,  h 


(50) 


*  — 2Lj.  ]  .  [0vh  sin(Bvh)  cos(Bvn)  -  Byn  cos(Bvh)  sin(Bvn) 
(x-5) 


♦  2  (2^)  cos(Bvh)  cos(Bvn)]  *  [S^h2^ 

S2h2 

-  3^5  1  [Syh  sin(Bvh)  cos(Bvn)  -  Svn  cos(Bvh)  sin(Bvn)  +  2  cos(Bvh)  cos(8vn)]> 


_ 2_ 

(x-0 


+  - 


2 

«n 

7' 


}  cos(ann)  sin(anz) 


Finally,  integrating  once  with  respect  to  x  and  z  one  finds* 


^  I  |  |  v(5,0,n>  Hjtlx-Cl ;n»z]  dndf- 

V”1  crack 
faces 


CD 

I 

|  I 

{±  v(g,o,n))  H2 [ | x-C I ;n*z)  dnde 

n=l 

crack 

faces 

1 

\ 

{+  v(5,0,n)}(2j-p)  dr|dS 

crack 

faces 

(52) 


-  -  (jq)  xz  ;  |x|  <  c  ,  |z|  <  h  . 

We  have  reduced,  therefore,  the  problem  to  that  of  the  solution  of  a 
two-dimensional  singular  integral  equation  for  the  unknown  function 
v(£,0,n).  This  solution  will  be  discussed  in  a  subsequent  paper. 

It  is  interesting  to  note  that  equation  (49)  is  also  applicable  to 
planar  cracks  of  arbitrary  shape  that  lie  on  the  x-z-plane  and  are 
symnetric  with  respect  to  both  x  and  z  -axes**. 

*The  reader  should  notice  that  the  function  v(5,0,n)  has  a  ~  sign  also. 
**The  same  method  of  solution  may  also  be  used  in  order  to  derive  a  much 
more  general  integral  equation  which  applies  to  any  arbitrary  crack  shape 
or  void.  This  matter  is  currently  under  investigation  and  the  results  will 
be  reported  in  another  paper. 


Perhaps  it  is  instructive  to  point  out  some  of  the 
advantages  of  the  present  formulation  over  that  of  reference  [4],  These 
are: 

(i)  we  are  seeking  the  solution  of  one  integral  equation 

(ii)  the  unknown  function  is  real  and  has  physical  meaning 

(iii)  the  unknown  function  can  be  related  directly  to  experi¬ 
mental  observations 

(iv)  the  formulation  applies  to  a  large  class  of  planar 
crack  problems 


V 

h 


i 
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6.  Solution  to  Navier 1 s  equations. 

Without  going  into  the  mathematical  details,  we  may  now  write 
the  displacement  functions  u^  ,  and  w^  in  terms  of  the 

unknown  function  v  (C,0,n)  ,  for  |x|  <  c  : 


u<c>  -  ± 


i  T  fh  vC5,o,n){  x~£  7 — T  -  sfl  i-  [lyl2  _ _ _] 

^i_cJ-h  ^  (x-5)z+|y|z  lfx^n2+lvlz 


Cx-O  *|yl‘ 


.  1_  ,2h2  .2  _2,  32  rx-K  , 

3E  CT“  z  n  J  — 7  [7 "I'  ,2*  ^ 

3x  (x-5)  +|y|A 

1  f c  fh  a  a 

±  ™  L  i-h V(C,M)  *  £  w  *  * 


(53) 


,C<0 


4  ^  C  l.h  T(5’M)  •5Jdn  d; 


1  r2h2  .2  2,  a2 


+  3m  ^T~  “  2*  "  "^7  2TT7^  41  ^ 

5  3xz  (x-e)%|yjz 


■  ®  l!c  Lh  TK’°>n)  •  st  cs^ 41  ** 


(M) 


ST  fc  ^.M)  $  *  « 


J 


*■1 
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wCc)  ■ 1 4  l!h  CV(E-M)  •  * .  &  CfciTF1 115  dn 

±  ^  |  1°  v(c,0,n)  ^  (M-M)  d5  dn  , 


(55) 


where  for  sinplicity  we  have  adopted  the  following  definitions: 


-  yvW.|yl2J 

N  =  2  J  - 2 — ,, - cosa^z  cosa_n 


(56) 


n-1 


n  n 


M  5  ^  v-1  1  ^  *  C°s2®vh)  COsSvh  cosV 


cosSyh  sin0vz>  {  -  i-jr  X^/fr-O^Iyl Z]  • 
ix 


I6vh  sin3vh  cosBvn  -  Byn  cos6vh  sinfivn+  2  (^)  cos0yh  cos0yn] 


(57) 


+  ev  K0(s/(x.0MxVi  (  Syh  sinSyh  cos3yn  -  Byn  cos3yh  sin3yn  ♦ 


+  2  cos3yh  cos0vn)) 


and 


.  # 


fc. 


1 
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«  2  .  ^  ^ 

I  *T“  cosB  h  cosg  z  {  -  -iy  K  [$  /(x-02+|y(Z] 
v=l  e^>  v  v  ax 

•  [  Bvh  sin0vh  cosS^n  -  Byn  cosByh  sinByii  +  2(~)  cosByh  cosg^n  ] 

_  (58) 

+  ev  V^/(x'^+M  ‘  f  Kh  sin0\,h  cosfy>  ‘  ®vn  cosByh  sin8vn 

♦  2  cosByh  cosBvn]}  • 

In  view  o£  the  above,  it  appears  that  the  solution  may  not  be 
separable  either  in  cylindrical  or  spherical  coordinates. 

Finally,  one  may  express  the  total  strain  energy  stored  in  the 
system  to  be: 

h  c 

W  -  -  i  f  [  {(v+  -  v")o  }  dxdz  .  (59) 

2  J.J.,  yy-o 


t 


3 


7.  Discussion 

Although  we  have  put  forth  a  considerable  amount  of  effort  to  solve 
the  double  singular  Integral  equation,  we  have  not  as  yet  been  successful 
In  recovering,  explicitly,  the  unknown  displacement  function  v  (x,  o,  z), 
valid  throughout  the  thickness  of  the  plate.  This  Is  an  extremely  difficult 
problem  where  physical  Intuition  can  be  misleading. 

At  present  we  have  developed  two  methods  which  In  principle  should 
give  us  the  desired  solution.  Unfortunately,  In  order  to  recover  the 
corner  singularity,  one  Is  forced  to  sum  up,  analytically,  a  double  series 
of  complex  el penf unctions.  This  Is  a  task  of  monumental  difficulty,  for 
the  algebra  Is  tedious  and  long. 
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APPENDIX  I. 

To  find  the  complete  homogeneous  solution  of  equation  (82)  of 
reference  [4],  we  proceed  as  follows. 

Assume  first  a  solution  of  the  form 

fwa*c)  -  ci-o2'2/m  gco  , 

where  G(?)  is  an  arbitrary  function  of  5  •  Next,  substitute  into 
homogeneous  difference-differential  equation  to  find 

(i-03"2/m  G'(0  *  (i+;)3'2/m  G'C-0  -  o  , 

from  which  one  may  now  deduce  that: 

■  Jo  f  U*;)3'2/m  tt)2"*1  .  c0 

or 

»  2n+2  - 

CCO  ■  lo  2F1  (-3*1  ,2n+2;2n*3;-c) 
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APPENDIX  II 

The  roots  of  the  equation  sin(2Bvh)  »  -  (28vh). 

The  equation  has  an  infinite  number  of  complex  roots  which  appear 
in  groups  of  four.  However,  as  it  was  pointed  out  in  the  text,  for 
this  analysis  only  the  roots  with  positive  real  parts  are  pertinent 
and  furthermore,  the  only  real  root  Bv  “  0  must  be  discarded.  Thus, 
if  we  define  the  roots  ®2,^4,^6> * * *  t0  the  complex  conjugates  of 
the  roots  3^,63,85,...,  then  by  setting 

28vh  *  \  *  %  v  *  1,3,5... 
and  using  a  Newtom-Rampson  numerical  method  one  finds 


B 

Xv 

1 

4.21239 

2.25073 

3 

10.71254 

3.10315 

5 

17.07337 

3.55109 

23.39836 

etc. 

3.85881 

Furthermore,  the  asymptotic  behavior  of  the  roots  for  large  v,  i.e., 
for  v  *  15,17,19,...,  is  given  by  the  following  simple  relations 


xv  -  (v  +  J-)ir 

yv  -  cos  h_1[(v  ♦  . 


_ 


Ob 


III!! 

'V'.  — . . . 


■I 

8888 

SI 

2h 


Figure  1.  Geometrical  representation  of  an  infinite  cracked  plate 
with  thickness  2h  and  crack  length  2c  . 
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PARTIAL  THROUGH  CRACK  IN  A  PLATE 
OF  FINITE  THICKNEXX 
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V 
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(P) 
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U  . 

v  , 

w 

x»  y, 

a 

n 

z 

V  rv*  Sn 


Yc  jng's  modulus 
E 

2(1+V) 

thickness  of  the  plate 

1 

v 


Fourier  transofrm  parameter 
displacement  functions 

displacement  functions  due  to  the  comple¬ 
mentary  problem 

displacement  functions  due  to  the  particu¬ 
lar  problems 

rectangular  cortesian  coordinates 
nir/h  n  -  1,  2,  3,  •** 
roots  of  the  eq.  sin  (f?vh)  ■  (g^h) 
roots  of  the  eq.  sin  (y^h)  »  -(y^h) 
coefficients  as  defined  in  text 


e 


3u  3v  .  3w 
3x  3y  3z 


v 


Poisson's  ratio 


*y’ 

Txz’  T: 

yz 

* 

r<c\ 

a'c>. 

a(C>, 

a<c> 

z  * 

xy  ’ 

xz  * 

yz 

rfP), 

c<P>, 

a(P), 

a<P> 

2  * 

xy  * 

xz  ’ 

yz 

stress  components 


*■  stress  components  due  to  the  comple¬ 
mentary  solution 

■  stress  components  due  to  the  particular 
solution 


uniform  applied  stress 
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I.  INTRODUCTION 

One  of  the  problems  In  fracture  mechanics  which  apparently  has  not 
received  extensive  theoretical  treatment  is  that  concerning  the  effect  of 
a  partial  through  crack  upon  the  stress  distribution  in  a  plate  of  finite 
thickness.  This  lack  of  interest  is  primarily  due  to  the  fact  that  three 
dimensional  problems  present  mathematical  complexities  which  are  substan¬ 
tially  greater  than  those  associated  with  plane  stress  or  plane  strain. 
However,  it  is  now  possible^ to  study  this  complex  phenomenon  which  has 
defied  researchers  for  some  time. 


II.  FORMULATION  OF  THE  PROBLEM 

Consider  the  equilibrium  of  a  homogeneous,  isotropic,  elastic  plate 
which  occupies  the  space  |x|<°°  ,  |y|  <  ®  ,  0  <  z  <  h  and  contains 

a  plane  crack  in  the  xz-plane.  The  crack  is  elliptical  in  shape  and  is 
defined  by  the  inequality 

<7>2  +  (7)2<  1  <« 

c  a 

The  plate  faces  z«0  and  z«h  are  free  of  stress  and  constraint.  Load¬ 
ing  is  applied  by  the  periphery  of  the  plate  |x|  ,  |y|  -►  00  and  is  given 
by 


ax  "  Txy  "  Tyz  "  0  *  ay  "  ao  • 

In  the  absence  of  body  forces,  the  coupled  differential  equations 
governing  the  displacement  functions  u  ,  v  and  w  are: 


(2) 

(3) 

(4) 


*See  references  [1,2  ] 


where 


a  _  9u  .  3v  .  3w 

0  “3J  +  W  +  ^ 

and  the  stress-displacement  relations  are  given  by  Hook's  lav  as: 

-  2G  { IS  +  J^2  J 

oy-20fg  +  £> 


3w 


-  2G  f  IS  +  S=2  J 


v0(g+g> 


As  to  boundary  conditions,  one  must  require  that  at: 


z  * 

0  : 

Txz 

’  Tyz  "  az  "  0 

z  ■ 

h  : 

Txz 

“  Tyz  "  °z  -  o 

(J)2  + 

<i>2 

<  1 

»  y  * 

0~  :  T  ■  t 
xy  yz 

|y| 

*+•  00 

and 

all 

x  :  Txy  "  Tyz 

1*1 

■+■  00 

:  o  »  T  ■ 
x  xy 

m  aym  0 


zx 


(5) 

(6) 

(7) 

(8) 
(9) 

(10) 

(ID 

(12) 

(13) 

(14) 

(15) 

(16) 


It  is  found  convenient  to  seek  the  solution  to  the  crack  plate  prob¬ 
lem  in  the  form 

u  -  u(P>  +  U(C)  etc.,  (17) 

where  the  first  component  represents  the  usual  "undisturbed"  or  "particular" 
solution  of  a  plate  without  the  presence  of  a  crack.  Such  a  particular 
solution  can  be  easily  constructed  and  for  the  particular  problem  at  hand  is 


* 


r«  -- 


[1-  (m-1)2] 


2,  qp(m~?> 


-  _  (*-2)2  (^)  z  (20) 

A  -  (a-l)s  -  3  (m-1)  +  2  (21) 

MATHEMATICAL  STATEMENT  OF  THE  COMPLEMENTARY  PROBLEM 

In  view  of  Che  particular  solution,  we  need  to  find  three  functions 
(C)  (C)  (C) 

u^  (x,y,z)  ,  v'  '(x,y,z)  and  w  (x,y,z)  ,  such  that  they  satisfy  simul¬ 
taneously  the  partial  differential  equations  (2)  -  (4)  and  the  following 
boundary  conditions: 

at  (f>  +  (f)2  <  1  ’  M  "  0:Txy(C)  “  Ty*(C>  "  °  ’  ay<C>  (22) 


at  z  *  0 
at  z  *  h 

at  V  x2  +y2  -*■  00 
METHOD  OF  SOLUTION 


”,/c>  -  V<c>  -  °,<c>  -  •  ») 

=  T«<C>  -  V2«>  -  ..«>  -  0  (24) 

:  u(C>  ,  v(C)  and  w^  are  to  be  bounded.  (25) 


In  constructing  a  solution  to  the  system  (2)  -  (4)  we  use  the  method 
described  in  reference  [l]  to  recover  the  following  ordinary  differential 
equation  of  the  Independent  variable  z  : 

*  »2  +*  3!>u<c)  +  V.*'®  +  3.>S<C>  ■ 0  <26> 

+  (I>2+  ^32>’'<c>+<^2  3.3i>“!c>  +  <^2  3*>  If—  • 0  <27> 


,2m- 2.  d2w(C) 

(T3T)  IP — 


♦<*  V  ^  »,>  ir1 *  °2»(c>  - 0 


‘  m-2  i'  dz 

a  t\2 


where  the  symbols  of  differentiation  3j  ,  32  ,  Dz  are  to  be  interpreted  as 
numbers. 


..  > 
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Upon  integrating  the  above  system  subject  to  the  initial  conditions*. 


(C) 


uo  *  v 


(C) 


(C) 


.CO 


(C)  du""  I  dvv 

,  WN  -  W.  ,  - -  “  U4  ,  J -  » 

»  0  *  dz  m  •  Am 


i  9  dz 


Cc) 


—  ”  w1  ,  for  z  *  0  , 
dz  o 

one  has  after  a  few  simple  calculations** 

“<c)  -  -  27^1,  ^  V.  *  “•  <-> 

“.‘a (.-2)  “ln.W  »,». 

’<C)  -  igife  ^  V.  -  u^T)  ■  - 


where 


+  cos(zD)  v _ S _ zstnSzDl  3  0 

-rcos^zu;  vfl  2(m-2)  D  2 


.(C) 


cos(zD)  w.  4- 


0  -  2(m~  Dz  Sin  (zD)  wo 


“A^e.+T (^2)  [^-Zcos(zD)]  ec 


0«  •  fri  (V„  +  32\> 


(29) 


(30) 


(31) 


(32) 


(33) 


Finally,  in  order  to  satisfy  the  boundary  conditions  (24  we  require  that 


[mhD  sin  (hD)31J  u0  +  [mhD92sin(hD)  J  v0 -mD  [  sin(hD)  -  hD  cos  (hD)]w#  ■  0  (34) 

[SinXhDl  (_  32_  (m  _  Dp2)  -  rah  cos  (hD)9*]  u0  +  9^ 

(35) 

-  mh  cos(hD)9j32]  vfl  +  [mltfj  Dsin(hD)]  w0  »  0 


*u0  ,  v0  ,  w0  ,  uj  ,  vj  ,  wj  are  arbitrary  functions  of  x  and  y 
**Hote  that  in  equations  (30)  _ (32)  we  have  let 

u  -  -  9  w 
o  l  o 


v  .  -3  w 

0  2  0 


in  order  to  simultaneously  satisfy  the  boundary  condition  (23) . 


-[mh  cos  (hD)3j32  +  sln^hDl  3j32]  u#  -  [mh  cos  (hD)  3?  +  g* 


+  (m-  1)  sin  (hD)  D]  vfl  +  [  hn  sin  (hD)  D  3,]  w.  -  0 


dn  d  1 2  d  1 3  uo  0 

d2 1  d22  d2  3  V0  *  0 

d 3 1  dS2  d  3  3  W0  0 


where  the  differential  operators  d^  are  defined  as 

d  ■  mh  3,  D  sin  (hD) 
n  1 

dJ2  ■  mh  32  D  sin  (hD) 

dj  3  -  -  mD  [sin(hD)  -  hD  cos  (hD)] 

d21  [sin(hD)  (3* -mD2) -mhD  32  cos(hD) ] 

d22  m~j)  f  sin(hD)  +mhD  cos  (hD)]  3j32 

dJ3  ■  mh3jD  sin(hD) 

dJ2  «  -  [sin(hD)  +mhD  cos  (hD)]  i  3132 

d32  “-J[sin(hD)(mD2-3j)+mh3*D  cos(hD)] 

dJ3  ■  mh  32  D  sin(hD)  . 

Keeping  in  mind  that  the  differential  operators  3j  ,  32  ,  D2  obey  the 
same  formal  rules  of  addition  and  multiplication  as  numbers,  the  solution  of 
system  (37)  is  given  by 

u0  “  Xi<*.y> 

v0  -  x2(*»y>  (3 

»0  ■  X|(*»y) 

where  the  unknown  displacement  functions  Xj  »  X2  »  X3  satisfy  the  differ¬ 


ential  relations 


i  *  1.2,3 
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dJl  d12  d13 


Q-  d21  d22  d23  -  m2(m-l)D*sin(hD)  (h2D2  -  sin2(hD)} 


d3  1  d  3 2  d 3  3 


We  construct  next  the  following  integral  representations  for  u.  ,  v 

9  0 

and  wQ  which  have  the  proper  behavior  at  Infinity 

»,(«,**>  •_sW  + 1  4’> « 

0  \T»  (42) 

+  l  rO>  +  l  sCOe-^W'ly!  }sin(xs)ds 

v=i  ^  n-i  n 

v0(x,y±)  -  +r(CP2+  |t|Q2)  e“8|y|  +  l  e“’^Tlyl +J  &<2> 

Vmi  V_1  (43) 

+  J  S^2^  e  _y/®  "^n^ ^  }cos(xs)ds 
n- 1  n 

».cr,^)  -r«p,+iyi«,).-|y|  +  I  i 

"0  3  3  VI  V  VI  V 

(44) 

+  £  S<3>e  2+ctn  ^  }cos (xs) ds 

n-i  n 

The  ±  signs  refer  to  y  >  0  and  y  <  0  respectively,  (a*  1,2,3, •  ••) 


and  ,  Yv  are  the  roots  of  the  equations 


sin(0vh)  -  (Bvh) 


sin  (Yyh)  ■  ”  (Yvh)  .  (46 

The  equations  have  an  infinite  number  of  complex  roots  which  appear  in  groups 
of  four,  one  in  each  quadrant  of  the  complex  plane  and  only  two  of  each 
group  of  four  roots  are  relevant  to  the  present  work.  These  are  chosen  to 
be  the  complex  conjugate  pairs  with  positive  real  parts.  The  only  real 
roots  ■  Yv  ■  0  must  be  ignored. 


- _ V 
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Finally,  an  examination  of  the  solution  shows  that  the  unknown 
functions  F  ,  Q  ,  R  etc.  are  not  all  independent.  Assuming,  there¬ 
fore,  that  one  can  differentiate  under  the  integral  sign  and  inserting 
equations  (42)- (44)  into  (37)  one  finds 


Q2  -  -Qx 

(47) 

(1  +  m)  a  (P1+P2)  +  (3m-l)  Qx+2smh  Q,  -  0 
d-cos(3vb)) 

sRv 

(48) 

(49) 

hR^2>  -  (l-co<8vh))Rv 

(50) 

V°  "  RV 

(51) 

<1  +  cos  V  Rv 
'v 

(52) 

hR^2)  -  "'yT  ■  (l  +  cosyvh)  S:v 

(53) 

-  * 

(54) 

V  V 

-  0 

(55) 

n 

,  > 

s(i)  - - a  s 

(56) 

n  s  n 

s(2>  -  s 
n  n 

(57) 

In  order  to  facilitate  our  subsequent  discussion  it  is  found  convenient  at 
this  stage  to  summarize  our  results: 

(i)  complementary  displacements*: 


“<c)  iV(pi  +*•>+ mv* 


-sfyl 


*It  can  be  shown,  that  in  order  to  satisfy  the  remaining  boundary  conditions 
must  vanish.  This  information  is  used  when  writing  the  complementary 

displacements  and  stresses. 


8 


•  1  -  cos8  h 

'v-i  8  ^  f ' r"  fife'  J  f co8gvz  - loSl)  gvz  8ln  Kz I 


“  2(^-1)  ^  [«dn(Pv2H^2  gvzc08®vz)1} 

cos(Yvh))‘ 

>^5 —  1008  M“  275^7  V 8lnV>l 


"  (l  +  co»{y  h})' 

'  8V‘ 

v*l 


(58) 


a- 2 
2(m-l) 


y-  t  sintf^  +^2  V  CO80rv4 ]> 


.(C) 


‘  Jx  SnC08<V)8"’/^T^ly|  }  •!»(*•)«»■ 

"  *J*o  {[P2'8zP3"2fe^T  z2sQ1-T(fc  *2*2CPX  +Pa)  -  lylQj] 


+  I  s  cos  fa  e  )e 


-/s2-Kxjy| 


n«i 


r  /t — r  ^1~cos(B  h)) 

+  fcoa(Bv25- evzsirf8vz)] 

P  [sin6vz  +~2  8vzcosB  z]}  e“,^'+®vlyl 


m-2 
2(m-l)  ft 


(59) 


"  a,  - 5-  (1+C08(Y  h)) 

\l  *-  ^  f-  V  '“•V-Icfe  V1"  VI 


nt-2 

2(m-l) 


p-  [sir<Yvz)+^2  Yv:c°8^]}e'^ily|}  cos(xs)d8  . 


.(C) 


„-8ly[ 


■Jf  {tP>  -si  *z  'p,+v -si ‘v 

•  (1-cos  (S  h)) 

+  j,  V - 57 - I2(SI7*1'<f,v‘)-2<SIT  V“*»v'» 

•  *  f  (l+cos(Yvh))  ^2 

+  v  {  Y^h  1 2(o-l)  *±nfrv^"2(m-l)  Yvz  C08<TVZ>] 

+COS  M+2(^TT  V  8inM  }e"/fl2+Yvlyl  |  coB(x8)dj  . 


(60) 
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6<c)-j;  <-2 <sr><,  t \ I cc.M 

-  6vhsin<Pvd]  «"’/®2+Svly!  (61 

•  &  n - 

+  ^f  Z  -h-[(1  +  C0*V®  co^rvZ)-Yvh8ln(Yv^]  e‘/82+Yvlyl}.coa(xs)d8  . 


(11)  complementary  stresses: 


t 

_Z£ 


(c) 


yZ  !•”  r  m  P  - 2  ,  ^ 

G  */,  {“S=I  /s  +e5  t— Th -  (»^+Bvzco*(8v^) 


-  S^z  sin(8vz)]  e~^a  +®vly  I 


(62) 


m  r  ^  rr~z  <1  +  coaM» 

~  o-l  Jj  Rv  /a  **%,  [ - 7h -  (Bll^v2)+Yvzcos(YvtD 


-  Yvz  sin(YvzJ]  e  ’^*YJyl  -j^  SQ  Oq  sin^e"^8^^  I  }cos(xs)ds  . 


I»  .r{J_  y  ,  (1"c°^ 

2G  J0  m-1  \r  g  h 


(sii^vz)+  Byz  cos(8v^)  -  Bvz  sin^zje-’^2*6^7 1 

•  v 

m  p  ^  (1  +  cosfr  ft)  /ITTi  1 

+  l  3  RV  I - Y  h  (sii^vz)+Yvz  cos(yv^D  -  Yvz  sin(Yvz»  e  8  +Yv'y' 


®  «42+a2 

+  1 

n-i 


as  s 

s  n  n  '  n 


lr(a_ z) e  ’^~^%lyl  }sin(xs)ds  . 


(63) 


(c) 


*V 


2G  "/,  ^  2 (m-1)  h  f  (*  -  CO{#v^)  6vz  sln(6vz)-  sln(6vz) 


+  8vh  8^z  cos(8vzJ]  e“  3  +®v^y^ 


*v 


(64) 


+  2C^U  Jml  T r  (1+co<V^  V  8inM"  Yvh  8in^Yvz) 

+  Yvh  V  C0<Yv^  8  ^  }cos(sx)ds  . 
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k~-J C  {f».- 

+  ~J  I  X  E  (1  -  co^Pv^- 6vh  •in(6vdl  e_,/*  +e,^y 

®  ra- co^b  ks> 

+  vix  Rv  (s2+8v)  [ - h -  f  C°8<M*  2(fe  6vzsil<0v4 

"  2^11  [  aia(W+^2  0v2  coaK$ 

+  f  S  cos  (a 

“  n  n  '  «  ' 


+  I  R 

V-l 


m-  2  _1_ 

2  (m-1)  8, 


n  n  □ 


*  X  rr~2 

+  T  [  (1  +  co^v^  c°^-Yvhsln(Yv4]  e“  s2+Yily| 

®  _  r  (l+co^  h) 

+  X  K  (s2+y^  — y^i  1  cosM- I(£it  v  slnM 


■  Y^  [  Sin^V')+i^2  V  cosM]  e"/8Z+Yvlyl|-cos(xs)ds  . 


T  lc) 

-ft  -+  J,“{[s(P1-P2)  +  2s 2 z  P3  +-2r-  z2s2Q1+2|y|sQ1  -Qjle'8!7! 


-  I  s*/^V 

\T1 


-  (1  -  cosfe  h& 

’v  V - jph - (  2  c08M'^1  0v*  sln^0v^ 


<sin(M+^  e„z  cos^.41  «“  8  +Bv'yl 


m-2  Mv' 


®  _ _  (l  +  co^Y.$) 

"  J, 8/8'‘^ K  [ — —  (2 co8M~^i  v 8iniyv*> 

-  El  ^  <  siniV>^  cos(Yvz>]e'/s2+Yv|yl 

®  (2s2  +a2)  ^5 — Ti  I 

-  I  - - — —  S  cos(a  z)  e"  8  ^  }  sln(xs)ds  . 

««■  i  *  ti  n 


f. 


-C{IsPi+s2zP3+5«:z282qi  +ysQi]  e"s|y| 

®  pci  -  cosffi  rt) 

-  X  a\  — Fh —  [  C0SM- 6vz  sinM 


m-  2  _1 

2  (m-1)  6 


j-  [  ffln(Svz)+  Svz  C08M]  C'^W 


+  ~j  I  [(1-  cos(6v^)  co^vz)-  6vh  sin(Bv^  e“  s  +gJyl 
00  „  Ri  +  cosCr.h) 

"X  [ — y$r~  tcos^V>T(fe  V  sinM 

-  T&lT  >T  [silM+^2  V  cosM  e-^X>yl 

+  Jj  T  [(1+  C0K\$>  coa(v^  Yvh sin(Yvz)]  e”  8 

-  y  /s2+ct2  S  co^a  z) e  8  ^n^  }cos(xs)ds. 

,  n  n  x  n  * 


i  i 

\  '4 


By  direct  substitution,  it  can  easily  be  ascertained  that  the  above 

complementary  displacements  satisfy  Navler's  equations  and  furthermore  the 

(c)  (c)  (c) 

corresponding  stresses  CP  ,  ,  T ^  do  vanish  at  the  plate  faces 

z  »  0  and  z  ■  h  . 

Moreover,  to  satisfy  the  continuity  conditions,  one  must  require 


I  8 
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c{ 


[p2  -  sz  P3  -  2(m_1)  z  sQ1  -  2 (m-1)  2  8  (P1  +  V ] 


+  7  s  cos (a  z) 

•  n  n 

n*l 


+  I  \  /s2+^2  {• 


(1-cos  (6  h)) 


v-1 


2,  [cos(f5vz)  -  2^iy  V  sin(Bvz)] 


6vh 


2 (m-1)  T  lsia(Kz)  +^2  $vz  co8^vz)l> 


(68) 


-  /2 — 5  (l+co8 (Y  h)) 

,Rv  /s  %  { - T -  [cosCYvZ)  -  2(^]J  V  sin^.z>] 

*  Y,*n 


"  2 (m-1)  T-  tsin<\z>  +^2  V  cos<V)]}}  cos(xs>ds 


,  ,  2  2 
+  tTc)  (x.o.z)  ;  ~  +  ~  <  1 
c  a 


2  2 

.  i_  +  >  1 

’2  2  ’ 
c  a 


.. 


]3 


c{ 


'*'!  +  ^  P3  +  275017  *V(!l  +  2(5017  *V(P1  +  P2> 


- 1 


1-0080,  h 


V»1 


- )  [cos(0vz)  -I7— -  6vz  sin(evz)] 


2(£l if  [sin(evz)  +  ^2  Bvz  cos^vz)]1 


(69) 


•  pj — t  ~  (l+cos(y  h)) 

-  ^  /s  +  3  Rv{ - -J- -  [cos  (Yvz)  -  yvz  8in(Yvz)] 


-  [sln(V>  +  ~  V  cos(V)]} 


2  2 
®  S  -KX 

I  — “  S_cos(a_z)^  sin(xs)ds 
n-1 


s  a  a 


-  3u(c)  2  2 

+  (fr-Vo  •  -T-1 


^  +  ■£»•  <  i 


c  a 


2  2 

2  2  X 

c  a 


and: 


I  t"-i  8  I  8 
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C  {[sP3  -  ==r  *2z  <pi +  p2>  -  “V 
/-= — s-  (1-cos (3  h))  , 

/shfi  ^ ^ -  [j^y  sin(0vz)  -  evz  cos(evz)] 


+  cos(6vz)  +jf~J  Bvz  sin(3vz)} 


(70) 


/-* — T~  (1+cos  (Y  h)) 

1's  *f-  K{ — ^ “ln(Yvz)  •  20irfT  V  “*(V>1 


+  "•  V>  +  2(^1)  V 


v»} 


Y  z  sln(y,  z)Jf  cos(xs)ds  ■ 


r-  3w(c\  x2  z2 

+<— Vo  ;  7  +  7<x 


2  2 
~  +  —  >  1 
c  a 


Thus,  by  Fourier  inversion,  one  has: 


■  t — »  8  ■  t-— i  8 
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!P2  '  *sP3  -  2^1)  *  2(fe  zV<Pl  +  VI 

CO 

+  I  Sn  coste^a) 

p — =•  (1-cos (Bh)) 

r^/s  +gy  { - - -  [cos(evz)  -  gvz  sln(3vz)] 

L  V1 

“  if^T)  IT  [sin(evz)  +  6vz  cos(Bvz)]} 

-  p — J  <1+cos(Y„h)) 

L  Vs  +Yv  { - ~r -  l*os(yvz)  -  2(SI3J  Yvz  Sin  Yvz)] 

'  2(£lT  ^  [sin(Yv2)  +  ^ 2  V  cos(V)]}  " 

§ 

-+|/0C  v(c)  (5,0,z)  cos(s5)dC 


(71) 


I 


'•pi +  *2*  P3 + ”s 


a(*>-0 


•***  0\)] 


®  /"5 — 5  1-cosB^h 

J.^  -%  8  Rv  « - )  tcos(Bvz)  -  2(^1) V  8in(Bvz)] 

20)  ^  tsln<8v*>  +  ^2  V  "*(V)1) 

*  /2 — 2~  -  (l+cos(Y  h)) 

J/a  %  8  V - ^ -  Ccoa(Yvz)  '  I^lTv2  8ln(V)] 

-  2^iT  ^  f8inCV)  +  A  V  C08(YVZ)1)  <72> 


y -  S  cos  (a  z)  - 

L.  a  n  n 

n-1 

,  S  a  (c) 

-  5  <!-?— ’Vo  8ln(8?)d5 

±4  £  <lr^Vo  “"<s5,d5 

+  /#c  v(c)  (5,0,z)  cos(sOd5 


(74) 
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from  which*, 

_  «  s 

8  S  -  +  ~  J“  /  C  V(5,0,n)  cos(aO  as  (a  n)dWn 
n  na  oo  n 

and 

.s0--2.  (P1  +  P2)Tii /•/%<“>  (C.o,«  • 

*  cos(s£)  d£d£. 


Finally,  we  would  like  to  solve  for  the  unknown  coefficients 
and  Rv  .  To  a  ccotnpllsh  this,  we  proceed  as  follows.  Using  eq. 
(71)  and,  upon  differentiating  w.r.t.  z  once,  one  has 


-  SP3  '  “Q1  *  ”2  <P1  +  V 


"  P) — J  (1-cos (Bh))  -  , 

+  v-1  { - V -  [  '  ifef)  S±n(8vZ>  "  2(£T)8v2 

+  [  -  cos(Bvz)  +  2  (to-lT^vz  8in<V)]} 


+  I  R  t 

M  *  M 


v-1 


(l+cos(Yvh))  3m_2 

"  1  ~  sin(V>  -  2G?I5V 


Yvh 


+  [  -  cos(Yvz)  +  2(5^T  V  sin<V)1} 

2  ,Sr  Sv(c) 

-  ±  cf)  /0C  (U—)y-0COs(a5)d5 


*  Notice  that  S  is  a  function  of  h  now. 
c 


(75a) 


(75b) 

*v 

(74)  into 


cos(&vz)] 

(76) 

cos(Yvz)] 


at; 
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Utilizing,  next,  the  orthogonality  condition*,  it  is  possible  to  determine 

3v(c) 

the  coefficients  and  in  terms  of  the  function  (3z  Vo  * 

Thus 


S  .  (c) 

%  -  C  l‘ 


y_0  B(Svh,6vn)  cos(s5)dC 


(77) 


and 


%  ■  *  « 


3v 


(c) 


k3n 


"Vo  B<Yvh,Yvn)  C08(85>d5 


(78) 


Finally,  inserting  into  equation  (65)  one  reduces  the  problem  to 
that  of  the  solution  of  a  double  singular  Integral  equation**,  i.e. 


// 


crack 

faces 


,3v 


(c) 


k30 


y-0 

aO 

'  2g 


H[x  “  (i^yiYy!2,Tl]d^dTl 


(79) 


where  the  kernel  H  consists  of  the  sum  of  three  infinite  series  of 
the  type  found  in  the  through-the-thickness  crack. 

The  explicit  solution  of  this  double  singular  integral  equation, 
will  determine  the  displacement  and  stress  fields.  Unfortunately, 
we  have  not  been  successful  in  extracting  the  solution  to  the  equation 
explicitly.  It  appears,  however,  that  the  solution  is  not  separable 
either  in  spherical  or  cylindrical  coordinates. 


*  See  reference  [2]. 

**  Eq.  (79)  may  be  Integrated  to  give  another  singular  integral  eq.  with 
(c) 

v  (C.O.n)  as  the  unknown. 
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PART  III 

UNIQUENESS  THEOREMS  FOR  DISPLACEMENT  FIELDS  WITH  LOCALLY 
FINITE  ENERGY  IN  LINEAR  ELASTOSTATICS 
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INTRODUCTION 

The  Classical  Theory  of  Linear  Elastostatlcs.  The  fundamental 
problem  of  linear  elastostatlcs  Is  to  determine  the  equilibrium  displace¬ 
ment  field  that  is  produced  In  an  elastic  body  of  known  shape  and 
composition  by  the  action  of  known  body  forces  and  surface  tractions  or 
displacements.  In  the  classical  formulation  of  the  theory  the  displace¬ 
ments  and  stresses  are  required  to  be  differentiable  and  satisfy  the 
differential  equations  of  equilibrium  in  the  interior  of  the  body  and  to 
be  continuous  and  satisfy  the  prescribed  surface  traction  or  displacement 
conditions  on  the  boundary.  This  boundary  value  problem  has  a  history 
that  begins  with  A.  L.  Cauchy's  discovery  of  the  equilibrium  equations 
in  1822;  see  reference  [18,  p.  8].  The  uniqueness  of  classical  solutions 
for  bounded  bodies  with  smooth  surfaces  was  proved  by  G.  Kirchhoff  in 
1859  [12].  General  existence  theorems  for  classical  solutions  were. first 
proved  during  the  period  1906-1908  by  integral  equation  methods.  The 
principal  contributors  were  I.  Fredholm  [6],  G.  Lauricella  [17], 

R.  Marcolongo  [19],  A.  Korn  [15,  16]  and  T.  Bogglo  [2,  3].  More  recently 
G.  Fichera  has  proved  the  existence  of  classical  solutions  in  bounded 
bodies  with  smooth  boundaries  by  the  methods  of  modern  functional  analy¬ 
sis  [4,  5].  Thus  the  theory  of  the  classical  boundary  value  problems  of 
linear  elastostatlcs  is  essentially  complete. 

The  Need  for  a  More  General  Theory.  Unfortunately  the  classical 
theory  described  above  provides  an  inadequate  foundation  for  the  analysis 
of  most  of  the  problems  studied  by  applied  scientists  in  their  applica¬ 
tions  of  linear  elastostatlcs.  Examination  of  any  of  the  numerous  books 
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on  theoretical  elasticity,  beginning  with  the  classical  treatise  of 
A.  E.  H.  Love  [18],  reveals  that  most  of  the  problems  treated  in  them 
involve  unbounded  bodies,  such  as  infinite  plates  or  bars,  and/or  bodies 

having  sharp  edges  or  corners.  Moreover,  the  stress  fields  are  known  to 

» 

have  singularities  at  re-entrant  edges  and  corners.  Examples  of  these 
difficulties  can  be  found  in  the  theory  of  cracks;  see  I.  N.  Sneddon  and 
M.  Lowengrub  [22] .  It  is  sometimes  argued  that  the  classical  theory  is 
a  sufficient  foundation  for  applications  because  real  bodies  are  always 
bounded  and  boundaries  with  sharp  edges  and  corners  can  be  approximated 
by  smooth  ones.  However,  although  this  procedure  simplifies  the  problems 
from  the  viewpoint  of  the  classical  theory,  it  makes  them  Inaccessible 
to  techniques  such  as  separation  of  variables  and  Integral  transform 
methods  that  are  used  by  applied  scientiests.  Thus  the  real  issue  is 
whether  a  mathematical  theory  can  be  devised  that  is  sufficiently  general 
to  provide  a  foundation  for  the  analysis  of  the  singular  problems  that 
are  actually  studied  by  applied  scientists.  The  purpose  of  this  paper 
is  to  provide  the  beginnings  of  such  a  theory  comprising  a  formulation  of 
the  elastostatic  boundary  value  problems  that  is  applicable  to  bodies  of 
arbitrary  shape  and  corresponding  uniqueness  theorems. 

Remarks  on  the  Formulation  of  Boundary  Value  Problems.  A 
"formulation"  of  a  boundary  value  problem  is  a  definition  of  the  class  of 
functions  in  which  solutions  are  to  be  sought.  The  classical  formulation 
of  the  elastostatic  boundary  value  problem  was  described  above.  Many 
other  formulations  are  possible.  For  example,  the  continuity  conditions 
may  be  replaced  at  some  or  all  boundary  points  by  boundedness  or  inte- 
grability  conditions,  the  equilibrium  equations  may  be  required  to  hold 
in  a  weak  sense,  etc.  In  principle,  any  formulation  is  acceptable  if 
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there  is  an  existence  theorem,  stating  that  there  is  at  least  one  solu¬ 
tion  in  the  class,  and  a  uniqueness  theorem,  stating  that  there  Is  at 
most  one  solution  in  the  class.  In  practice  the  choice  of  a  solution 
class  turns  on  technical  considerations.  The  proof  of  an  existence 
theorem  is  facilitated  by  choosing  a  large  solution  class  but  uniqueness 
is  lost  if  the  class  is  too  large.  The  proof  of  a  uniqueness  theorem  is 
facilitated  by  choosing  a  small  solution  class  but  existence  is  lost  if 
the  class  is  too  small.  For  example,  Kirchhoff's  theorem  on  the  unique¬ 
ness  of  classical  solutions  of  the  elastostatic  boundary  value  problem 
can  be  proved  for  bodies  having  re-entrant  sharp  edges  but  in  this  case 
no  classical  solution  exists. 

The  Role  of  Existence  and  Uniqueness  Theorems.  A  pure  existence 
theorem  for  a  boundary  value  problem  demonstrates  that  the  properties 
chosen  to  define  the  solution  class  are  not  contradictory;  i.e.,  there 
are  functions  with  these  properties.  In  the  presence  of  an  existence 
theorem  a  uniqueness  theorem  shows  that  the  defining  properties  of  the 
solution  class  characterize  the  solution  completely.  However,  a  unique¬ 
ness  theorem  can  be  even  more  valuable  when  no  general  existence  theorem 
is  known.  In  such  cases  it  may  still  be  possible  in  certain  instances, 
corresponding  to  special  choices  of  the  boundary  or  data,  to  construct  a 
solution  in  the  chosen  solution  class.  A  uniqueness  theorem  then  shows 
that  the  solution  is  the  correct  one.  An  interesting  example  of  this 
occurred  in  the  theory  of  the  diffraction  of  electromagnetic  waves  by  a 
perfectly  conducting  circular  disk.  In  1948  J.  Meixner  [20]  proved  a 
uniqueness  theorem  for  this  problem  and  used  it  to  show  that  a  solution 
that  had  been  published  in  1927  was  incorrect.  Of  course,  in  the  absence 
of  a  general  existence  theorem  it  is  desirable  to  prove  uniqueness  in  as 
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large  a  solution  class  as  possible  since  this  facilitates  application  of 
the  uniqueness  theorem  in  specific  Instances. 

The  Boundedness  Question  for  the  Displacement  Fields.  Linear 
elastostatlcs  is  an  approximation  that  is  valid  for  small  displacements. 

If  the  displacements  are  bounded  then  by  suitable  scaling  they  may  be 
made  arbitrarily  small.  Hence  it  is  natural  to  make  boundedness  of  the 
displacements  a  defining  property  of  the  solution  class.  Indeed,  this 
property  has  often  been  employed  in  constructing  solutions  of  particular 
problems.  It  has  also  been  used  by  J.  K.  Knowles  and  T.  A.  Pucik  [14] 
in  the  formulation  and  proof  of  a  general  uniqueness  theorem  for  plane 
crack  problems.  However,  it  is  shown  in  this  paper  that  uniqueness  holds 
in  the  larger  class  of  solutions  with  locally  finite  energy,  without 
boundedness  conditions.  This  result  shows  that  the  boundedness  hypothesis 
is  redundant  and  the  boundedness  property,  in  instances  where  it  holds, 
must  be  derivable  from  the  other  hypotheses. 

Displacement  Fields  with  Locally  Finite  Energy.  In  this  paper 
it  is  taken  as  a  fundamental  principle  that  equilibrium  displacement 
fields  in  elastic  bodies  must  have  finite  strain  energy  in  bounded  por¬ 
tions  of  the  bodies.  Such  displacement  fields  will  be  called  displace¬ 
ment  fields  with  locally  finite  energy  (or,  for  brevity,  fields  wLFE) . 

The  equilibrium  displacement  field  corresponding  to  prescribed  body 
forces  will  be  characterized  among  all  fields  wLFE,  by  the  principle  of 
virtual  work.  The  class  of  displacement  fields  that  obey  these  two 
principles  will  be  called  the  solutions  with  locally  finite  energy  (for 
brevity,  solutions  wLFE)  of  the  elastostatic  boundary  value  problems. 

The  principal  results  of  this  paper  are  uniqueness  theorems  for  this 
class  of  solutions.  In  particular,  the  uniqueness  of  solutions  wLFE  in 
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bounded  bodies  is  proved  without  additional  hypotheses  concerning  the 
boundary  or  the  displacement  field.  The  uniqueness  of  solutions  wLFE  in 
unbounded  bodies  is  proved  under  a  growth  restriction  on  the  behavior  of 
the  stress  or  displacement  fields  at  infinity.  Moreover,  it  is  shown  by 
examples  that  a  growth  restriction  is  necessary  for  uniqueness. 

The  remainder  of  the  paper  is  organized  as  follows.  The  class 
of  displacement  fields  wLFE  is  defined  in  §1.  §2  contains  the  definition 

of  the  class  of  solutions  wLFE  in  homogeneous  elastic  bodies  of  arbitrary 
shape,  subject  to  prescribed  surface  tractions,  prescribed  body  forces 
and  prescribed  displacements  or  stresses  at  infinity.  The  regularity 
properties  of  solutions  wLFE  are  also  discussed  in  this  section.  S3 
presents  the  uniqueness  theorems  for  solutions  wLFE  of  problems  with 
prescribed  surface  tractions.  In  §4  the  methods  and  results  of  §3  are 
extended  to  the  other  classical  boundary  value  problems  of  linear  elasto- 
statics  including  problems  with  prescribed  surface  displacements, 
problems  with  mixed  boundary  conditions,  problems  for  inhomogeneous 
elastic  bodies  and  n-dimensional  generalizations.  55  contains  a 
discussion  of  related  literature. 


1.  DISPLACEMENT  FIELDS  WITH  LOCALLY  FINITE  ENERGY 


A  fixed  system  of  Cartesian  coordinates  is  used  throughout  the 
paper  and  points  of  Euclidean  space  are  identified  with  their  coordinate 
triples  (xj,x2,x3)  -  x  £  R3.  With  this  convention  each  elastic  body  in 
space  is  associated  with  a  domain  (open  connected  set)  (2  C  R*  that 
describes  the  set  of  interior  points  of  the  body.  The  closure  and 
boundary  of  £2  are  denoted  by  £2  and  3£2  ■  £2  -  £2,  respectively.  The  nota¬ 
tion  of  Cartesian  tensor  analysis  [11]  is  used  to  describe  the  physical 
variables  associated  with  elastic  bodies.  In  particular,  tensors  of 
various  orders  are  denoted  by  subscripts  and  the  summation  convention  is 
used. 


The  fundamental  unknown  of  elastostatic  boundary  value  problems 
is  the  displacement  field.  It  is  denoted  below  by  u^  »  u^Cx) .  The 
notation  u 


3u^/3xj  is  used  for  the  covariant  derivative  of  u^.  The 


strain  tensor  field  e^j(u)  associated  with  u^  is  defined  by  the  differen¬ 
tial  operator 


(l.D  eij(u)  ”  2  (ui,j  +  Jj,i* 

It  is  assumed,  following  G.  Green  [7  and  18,  pp.  11-12  and  95-99],  that 
for  quasi-static  isothermal  small  deformations  of  an  elastic  body  there 
is  a  positive  definite  quadratic  function  of  e^, 

(1.2)  w  -  j  ciju  eij  %JL  » 

such  that  for  all  KC  12 
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(1.3) 


e,,(u)  e..(u)  dx 


,K  cijki  eijQu;  ekT 

is  the  strain  energy  of  the  displacement  field  in  the  set  K.  The 
positivity  assumption  means  that 


(1.4) 


Cijk£  eij 


e^  >  0  for  all  e^  -  e^  j*  0 


The  stress-strain  tensor  is  uniquely  determined  by  w  if  the  natural 

symmetries 


(1.5) 


cijk£  “  cj iki.  "  ck£ji 


are  assumed.  The  stress  tensor  field  a^(u)  associated  with  u,  is  given 
by  the  differential  operator 


(1.6) 


Oij(u)  -  C±jU 


The  positive  definiteness  of  w  implies  that  a '  »  c^ ^  e^  has  a  unique 

solution  e±i  -  YijU  a u  and  w  -  \  o±j  e^  -  £  o±j  O^.  In 

particular, 

d.?)  WK  -  |  aij  (“)  eij  (u)  **  "  I  }R  YljkZ  °ij  <“)  ou(u)  dx 

is  a  functional  of  CT^j(u)  alone.  A  body  is  homogeneous  if  and  only  if 
cijk£  is  constant  in  It:  Is  isotropic  if  and  only  if  [11,  18] 


(1-8) 


cijkl  “  X  6ij  6kJl  +  v(5ik  6jJl  +  5U 


where  X  and  y  are  scalars  such  that  y  >  0,  3X  +  2y  >  0.  The  results  in 
§2  and  S3  are  formulated  for  the  case  of  homogeneous  anisotropic  bodies. 
In  §4  it  is  shown  that  the  uniqueness  theorems  hold  for  the  more  general 
case  of  inhomogeneous  anisotropic  media  with  bounded  uniformly  positive 
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definite  stress-strain  tensor.  This  means  that  the  components 

are  Lebesgue  measurable  and  there  exist  positive  constants  cg  and  ct_>  c0 


such  that 


(1-9)  e0  e^  e^j  <  e^j  ^  1  ci  eij  eij  *or  811  x  e  fl 

and  all  e4^  «  e^. 

The  most  general  uniqueness  theorems  for  solutions  wLFE  will  be 
obtained  by  making  the  class  of  displacement  fields  wLFE  as  large  as 
possible  subject  to  the  LFE  condition.  Hence  it  is  natural  to  define  the 
energy  integrals  W^(u)  to  be  Lebesgue  integrals  and  to  Interpret  the 
differential  operators  e^  in  the  distribution-theoretic  sense.  It  can 
be  shown  that  this  choice  has  the  additional  advantage  that  the  set  of 
displacement  fields  wLFE  is  a  complete  space  in  the  sense  of  convergence 
in  energy  on  bounded  sets.  It  was  by  using  such  complete  function  spaces 
that  Fichera  proved  the  existence  of  solutions  of  the  elastostatic 
boundary  value  problems  in  bounded  domains. 

In  the  remainder  of  this  section  several  function  spaces  are 
defined  that  are  needed  for  the  formulation  and  proof  of  the  uniqueness 
theorems.  In  the  definitions  fl  C  Rs  denotes  an  arbitrary  domain. 

The  definitions  are  based  on  the  Lebesgue  space 

(1.10)  L2(fi)  ■  {u:  fl  -*■  R  |  u(x)  is  L-measurable,  u(x)2  dx  <  “»} 
and  the  associated  spaces 

L^°C(ft)  *  {u:  fl  -*■  R  |  u  6  Lj(K)  for  every  bounded 

(1.11) 


measurable  KC  Q) 
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(1.12)  L2nt(ft)  »  (u:  fi  +  R  |  u  e  L2(C)  for  every  compact  C  C  8} 
and 

t  .  L2°m(ft)  -  L2(ft)  n  (u  |  u(x)  is  equivalent  to  0  outside 

(1.13) 

a  bounded  set} 

i 

It  is  clear  that  L2om(ft)  C  L2  (ft)  C  L2°C(ft)  C  L2nt(ft).  Moreover, 

L^0™^)  -  L2(ft)  *  L2°C(ft)  if  and  only  if  0  is  bounded.  Note  that  the 
Qqc. 

condition  u  €  L2  (ft)  restricts  the  behavior  of  u  near  9ft  because  the 
sets  K  in  (1.11)  can  be  any  bounded  open  subsets  of  ft.  The  condition 
u  £  L2nt(ft)  is  weaker  because  it  does  not  restrict  the  behavior  of  u  near 
9ft.  All  of  the  function  spaces  used  below  are  spaces  of  tensor  fields 
on  ft  whose  components  lie  in  certain  linear  subspaces  of  L2nt(ft) . 

The  space  L2nt(ft)  may  be  interpreted  as  a  linear  subspace  of 
L.  Schwartz's  space  V' (ft)  of  all  distributions  on  ft  [21].  Thus  functions 
u  £  L2nt(ft)  have  derivatives  of  all  orders  in  V (ft)  and  if 

(1.14)  A-  l  A  9lQl/9x“l  9x22  9x“j 

0<Ja|<m 

(where  a  -  (o1,a2,aJ),  |a|  ■  ax  +  a2  +  os)  is  a  partial  differential 
operator  with  constant  coefficients  then  Au  £  O' (ft).  The  notation 
Au€  L2nt (ft)  (resp.  I.foc(n) -La  (ft)  ,  L2°m(ft)  ,  etc.)  will  be  interpreted 
to  mean  that  the  distribution  Au  is  in  the  subspace  L2nt(ft)  (resp. 
I*2°C(ft),  L2 (ft)  ,  L2om(fi),  etc.).  If  Aj ,  A2 , •  •  • ,  AQ  is  a  set  of  partial 
differential  operators  with  constant  coefficients  the  following  notation 
will  be  used. 

(1.15)  L2 (Al ,A2 , • • *,AQjft)  -  L2 (ft)  n  (u  |  AjU  £  L2 (ft) , j  -  l,2,***,n} 
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(1.16)  LfOC(A1,AJ,***,An;fi)  -  L2°C(ft)  n  (u  1  AjU  e  lJ°C(«)  ,J  -  l,2,***,n} 

(1.17)  Ljnt(A1,A2,---,An;fi)  -  I*nt(ft)  n  {u  |  AjU  €  L2nt(ft),j  -  1,2, •••,n> 

(1 .18)  L2on<A1  ,A2 ,  •  •  •  ,AQ;fi)  »  Lf  “(B)  O  L2  (A,  ,A2 ,  •  •  •  ,An;ft) 

In  particular,  if  {Ai.Aj,*** .A^}  »  {S^/Bx®1  3x22  3x°3  |  0  <  |a|  <  m) 
the  following  notation  will  be  used. 

(1.19)  L?(ft)  -  L2(Aj,A2,»-tAn;a) 

(1.20)  L“*£°C(ft)  -  Lf°C(Al,A2,*“,An;n) 

(1.21)  L2’lnt(ft)  -  lJnt(Al,A2,  —  ,Ato;Q) 

(1.22)  L2,COm(«)  -  Lj°“(Al,Aa,-*,An;n) 

0  QP 

Notations  such  as  u^^  €  L2  (ft),  e^  e  L2  (ft) ,  etc.  will  be  inter¬ 
preted  to  mean  that  each  component  of  the  tensor  field  is  In  the  Indicated 
space.  With  this  convention  the  classes  of  displacement  fields  wFE  (with 
finite  energy)  and  wLFE  may  be  defined  as  follows. 

Definition.  A  vector  field  u^  on  ft  is  said  to  be  a  displacement 
field  wFE  if  and  only  if  it  is  in  the  function  space 

(1.23)  E(fl)  -  {u  |  uA  €  L£oC(ft),  e±j(u)  6  L2(ft)> 

Similarly,  u^  is  said  to  be  a  displacement  field  vLFE  if  and  only  if  it 

is  in  the  function  space  • 

(1.24)  E£oc(ft)  -  {u  |  u±  6  L£oC(ft),  e^j (u)  €  lJ°C(B)} 

Note  that  E£oc(ft)  «  E(ft)  if  and  only  if  ft  is  bounded. 


The  terminology  used  In  Che  definition  is  justified  by  the  obser¬ 


vation  that  if  the  stress-strain  tensor  satisfies  (1.9)  then  e^ (u)  6  L2(8) 
implies  Cfy(u)  e  1*2(8)  and  hence  u  e  E(8)  implies 

(i.25)  »n  ■  T  |n  «i3(u)  ‘  * 

0q/»  0 rs  £ 

Similarly,  if  (1.9)  holds  then  e^  (u)  €  L2  (8)  implies  o^OO  e  1*2  (8) 

f/j  n 

and  hence  u  6  E  (8)  implies 


(1.26) 


W. 


K 


2  Ik 


cij(u) 


eij  (u) 


dx  < 


for  all  hounded  measurable  sets  K  c  8. 

Each  of  the  function  spaces  defined  above  is  a  complete  space 
with  respect  to  a  suitable  topology.  Several  examples  of  this  will  be 
indicated.  It  is  well  known  that  L2(8)  is  a  Hilbert  space  with  scalar 
product 


(1.27)  (u,v)  -  f  u(x)  v(x)  dx 

'8 

SjOc  + 

Similarly,  E(8)  and  E  (8)  are  Frechet  spaces  [28]  with  respect  to  the 
families  of  semi-norms  defined  by 


(1.28) 


PK,E(U) 


(1, 


uA(x) 


ut(x) 


dx  + 


I* 


aij(u) 


Vu) 


*1 


1/2 


and 

(1.29)  Pk?e<u)  “  (|  “*(*)  +aij(u>  0ij^u^  ta) 

respectively,  where  K  is  any  bounded  measurable  subset  of  8.  In  parti- 

0  Qc 

cular,  if  8  is  bounded  then  E  (8)  *  E(8)  is  a  Hilbert  space.  These 
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completeness  results  play  no  role  in  the  uniqueness  theorems  given  below. 
However,  they  are  essential  for  the  validity  of  existence  theorems  for 
solutions  wLFE.  This  is  evident  from  the  proofs  of  Fichera's  existence 
theorems  for  bounded  bodies. 

j&oc 

In  Che  definition  of  E  (ft)  the  operators  e^(u)  defined  by 
(1.1)  are  interpreted  in  the  distribution-theoretic  sense.  Hence  the 

£o  Q 

condition  u  €  E  (ft)  does  not  necessarily  imply  that  the  individual 
derivatives  u^  6  Lj0C(ft).  However,  it  is  known  that  if  u  €  E£oc(ft) 
then  u^  €  L2(C)  for  every  compact  set  C  C  ft.  This  is  a  consequence  of 
Korn's  Inequality  in  the  form 

I  u 

Oq  /* 

which  is  valid  for  all  u  6  E  (ft) ,  all  bounded  open  sets  K  C  ft  and  all 
compact  sets  C  C  K  with  a  constant  y  •  y(C,K).  This  result  can  be 
derived  from  the  version  of  Korn's  inequality  due  to  J.  Gobert  [8). 
Moreover  if  ft  has  the  cone  property  [1,  9]  then  one  may  take  C  ■  K  in 

(1.30) .  Hence  in  this  case 

(1.31)  u  6  E^°c(ft)  -  u±  e  L2’£oc(ft) 

In  particular,  for  domains  that  are  bounded  and  have  the  cone  property 

(1.32)  u  €  E(ft)  •  u±  €  L|  (ft) 


i' L, 


(K) 


3 

l 

i.j-1 


(K) 


(1.30) 


“i.J *MC)  -  Y 


i-1 


2.  EQUILIBRIUM  PROBLEMS  WITH  PRESCRIBED  SURFACE  TRACTIONS 

In  this  section  elastostatic  equilibrium  problems  are  formulated, 
and  regularity  properties  of  the  solutions  are  discussed,  for  homogeneous 
anisotropic  elastic  bodies  of  arbitrary  shape  that  are  subject  to 
prescribed  body  forces,  prescribed  surface  tractions  and,  in  the  case  of 
unbounded  bodies,  prescribed  displacements  or  stresses  at  infinity.  The 
cases  of  prescribed  body  forces  F^,  zero  surface  tractions  and  zero  dis- 
placements  or  stresses  at  infinity  are  discussed  first. 

The  Principle  of  Virtual  Work.  Let  £2  C  R*  be  an  arbitrary  domain 

J?  oc 

and  let  u  €  E  (£2)  be  the  equilibrium  displacement  field  vLFE  correspond¬ 
ing  to  body  forces  Fj,  6  Lj0m(£2)  and  zero  surface  tractions.  Imagine  that 
the  equilibrium  is  disturbed  slightly  by  changing  u.^  to  uA  +  where 
is  a  field  wFE  from  the  set 

(2.1)  ECOm(£2)  -  E(£2)  fl  {v  |  e^Cv)  €  lJ°“<»)> 

Let  K  C  £2  be  a  bounded  measurable  set  such  that  e^(v)  is  equivalent  to 
zero  in  £2  -  K.  Then  V?K(a(u) )  and  WK(o(u  +  v))  are  the  strain  energies  in 
K  before  and  after  the  disturbance.  Hence  the  work  done  against  Internal 
forces  during  the  disturbance  is  Wj,(o(u  +  v))  -  Wj^(a(u)).  The  energy 
norm  of  v  can  be  made  arbitrarily  small.  If  this  is  done  and  terms 
quadratic  in  v  are  dropped,  in  keeping  with  the  linear  theory,  the 
difference  becomes 

(2.2)  o.,(u)  e,  (v)  dx  -  Work  done  against  internal  forces 
'£2  *1  1-1 

Moreover,  if  the  body  forces  are  constant  during  the  displacement  then 
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(2.3)  -  F,  v.  dx  ■  Work  done  against  body  forces 

No  further  work  is  done  during  the  disturbance  if  the  surface  tractions 
are  zero.  The  principle  of  virtual  work  states  that  the  true  equilibrium 
field  u^(x)  is  characterized  by  the  property  that  the  total  work  done 
against  the  Internal  and  external  forces  in  any  (small)  disturbance  of 
u^  consistent  with  the  constraints  is  zero  [23].  Thus  in  the  present 
case 

(2.4)  |  °ij(u)  eij^  ^  F £  v±  dx  -  ° 

for  all  v  €  Ecom(fl) .  This  motivates  the  following 

Definition.  A  displacement  field  u^  is  said  to  be  a  solution 
wLFE  of  the  equilibrium  problem  for  the  domain  £2  with  body  forces 
1?^  e  L j°m(fl)  and  zero  surface  tractions  if  and  only  if  u  6  E^OC(I2)  and 

(2.4)  holds  for  all  v  €  ECOffl(n). 

Necessary  Conditions  for  the  Solvability  of  Problems  with  Zero 
Surface  Tractions.  The  fields 

(2.5)  v±(x)  -  a±  +  bj  xfc  ,  x  6  R* 

where  aA  and  b^  are  constant  vectors  and  is  the  alternating  tensor 

[11]  satisfy  ey(v)  ■  0  in  Rs  and  hence  v  6  Econ(Rs).  In  particular, 
v  €  Ecom(£2)  for  every  domain  £).  It  follows  from  (2.4)  with  this  choice 
of  v  that  necessary  conditions  for  the  existence  of  a  solution  wLFE  are 

(2.6)  [  F.  dx  -  0 

(2.?)  j  (F±  Xj  -  Fj  xA)  dx  -  0 
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Physically,  these  conditions  mean  that  the  body  forces  F^  exert  no  net 
resultant  or  moment  on  the  body.  They  are  assumed  to  be  satisfied  in 
the  remainder  of  the  discussion  of  problems  with  zero  surface  tractions. 

Non-uniqueness  of  the  Displacements  for  Problems  with  Zero 
Surface  Tractions.  Equations  (2.5)  define  a  displacement  field  that 
describes  a  rigid  body  displacement  [11].  Moreover,  since  e^ (v)  ■  0  in 
R3  the  fields  (2.5)  may  be  added  to  any  solution  u  of  (2.4).  Physically, 
this  means  that  the  equilibrium  displacement  fields  are  determined  only 
up  to  rigid  body  displacements.  Hence,  the  natural  uniqueness  theorem 
for  problems  with  zero  surface  tractions  asserts  that  the  stress  and 
strain  fields  are  unique  while  the  displacement  fields  are  unique  modulo 
fields  of  the  form  (2.5). 

Bounded  Bodies  and  Displacement  Fields  wFE.  If  ft  is  bounded  then 
E  (ft)  *  E(ft)  and  every  solution  wLFE  actually  has  finite  total  strain 
energy  in  ft.  More  generally,  if  u  is  a  solution  wLFE  for  an  arbitrary 
domain  ft  and  if  u  €  E(ft)  then  u  is  said  to  be  a  solution  wFE.  The 
uniqueness  of  solutions  wFE  is  proved  in  §3  without  additional  hypotheses 
concerning  ft  or  the  displacement  field. 

Unbounded  Bodies  and  Equilibrium  States  with  Prescribed  Stresses 
or  Displacements  at  Infinity.  If  ft  is  unbounded  then,  in  general,  solu¬ 
tions  wLFE  in  ft  are  not  unique.  Simple  examples  of  non-uniqueness  are 
available  for  the  case  ft  ■  R3.  The  field  u^(x)  »  b^x^  with  constant 
b^  ■  bj^  4  0  is  a  solution  wLFE  in  R*  with  F^(x)  =  0  and  a^fu) 

*  cijkZ  bkA  *  0  since  eij(u)  "  bij  and  0ij(u)  eij(u)  “  Cijki  bij  bk 1  >  °* 
A  second  example  is  provided  by  the  homogeneous  Isotropic  plate  with 

domain  ft  ■  {x  |  Xj,x2  6  R, |x  }  <  h)  and  stress-strain  tensor  (1.8).  In 

this  case  uj  •  (X  +  2p)xlt  u4  -  (X  +  2y)x2,  u,  -  -2Ax,  defines  a 


I 
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displacement  field  in  ft  with  F^(x)  =  0,  zero  surface  tractions  and 
constant  non-zero  stress  field  0j l  ■  o22  ■  6Xy  +  4y2,  all  other  «  0. 
These  examples  show  that  uniqueness  theorems  for  solutions  wLFE  in 
unbounded  domains  cannot  hold  without  some  groxjth  restrictions  at 
infinity  on  u^  or  0^j . 

The  problem  of  finding  suitable  growth  restrictions  on  ui  or 
that  guarantee  the  uniqueness  of  solutions  wLFE  is  a  special  case  of  the 
classical  problem  of  elastostatics  of  finding  equilibrium  displacement 
fields  that  have  prescribed  stresses  or  displacements  at  infinity.  Many 
problems  of  this  type  are  discussed  in  the  treatise  of  Love  [18] .  To 
formulate  the  problem  with  prescribed  stresses  at  infinity  let 


(2.8) 


-  ft  n  {x  |  | x I  >  R} 


<30 

and  let  0^  (x)  be  a  stress  field  that  is  defined  in  S,,.  oo,  for  some  R, 
and  has  the  desired  behavior  at  infinity.  A  solution  wLFE  in  ft  is  sought 

OO 

such  that  o^jCuJCx)  is  close  to  0^(x)  at  infinity,  in  a  suitable  sense. 

OO 

One  possibility  is  to  require  that  O^Cu)  -  6  L^i^  a)  or> 

equivalently, 


(2.9) 


This  suggests  the 


(a(u)  -  a")  <  » 

R,w 


Definition.  A  solution  wLFE  of  the  equilibrium  problem  for  an 

OO 

unbounded  domain  ft  is  said  to  have  prescribed  stresses  0^  at  infinity 
if  and  only  if  (2.9)  holds  for  some  R  >  0. 

OO 

Solutions  wLFE  with  stresses  0^  ■  0  at  infinity  are  just  the 
solutions  wFE  defined  above.  Condition  (2.9)  is  correct  in  this  case, 
at  least  for  exterior  domains  where  the  stresses  generated  by  body  forces 
6  L‘on(ft)  are  known  to  satisfy  ©^(uKx)  -  0(|xj”2),  |x|  -*■  00  [13]. 
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To  formulate  the  problem  with  prescribed  displacements  at  infinity 
let  u*(x)  be  a  displacement  field  that  is  defined  in  B  for  some  R,  and 
has  the  desired  behavior  at  infinity.  A  solution  wLFE  in  ft  is  sought 

OQ 

such  that  u^x)  is  close  to  u^(x)  at  infinity,  in  a  suitable  sense.  One 

00 

might  try  the  condition  u^  -  u^  e  L^ft^  ,  in  analogy  with  (2.9). 

However,  this  condition  is  too  strong.  In  fact,  it  is  known  that  if 

CO 

u^  •  0  and  ft  is  an  exterior  domain  then  the  displacements  generated  by 
body  forces  F^  €  Lj0in(ft)  have  the  exact  order  u^(x)  ■  0(|x|  1),  |x|  ® 

[10].  Thus  a  weaker  conditi.n  consistent  with  this  estimate  is  needed. 

In  what  follows  the  condition 

(2.10)  lu  -  uV  .  -  £)( r),  r 

r»o 

is  used  where 

(2.11)  ,u,r  5  ”  jft  ui<x>  ui(x>  d*  » 

r ,  5 

(2.12)  B  ,  -  n  n  (x  I  r  <  |x|  <  r  +  6} 

r,o  -  ~ 

and  £  >  0  is  a  constant. 

Definition.  A  solution  wLFE  of  the  equilibrium  problem  for  an 

00 

unbounded  domain  ft  is  said  to  have  prescribed  displacements  u^  at 
infinity  if  and  only  if  (2.10)  holds  for  some  5  >  0. 

A  sufficient  condition  for  (2.10)  to  hold  with  u*  «  0  is 

(2.13)  ut(x)  -  0(jxr1/3)  ,  |x|  ® 

Of  course,  the  precise  order  condition  on  u^  that  is  sufficient  to  guar¬ 
antee  (2.10)  in  particular  cases  will  depend  on  the  geometry  of  ft  near 

Infinity.  For  example,  if  ft  »  {x  |  |x,'|  <  h}  then  fn 

*  “r,6 


dx  -  0(r) 
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r  -*•  «®,  and  u^(x)  ■  0(1)  is  a  sufficient  condition  for  (2.10)  with  »  0. 

If  fl  ■  (x  |  (Xj,x2)  €  G,  x3  6  R}  where  GCRJ  is  bounded  then 

/ft  dx  -  0(1)  and  u .  (x)  ■  (  | x J 1/2 ) ,  |x|  •*■•*>,  is  sufficient. 

r,o  1 

Ellipticity  of  the  Cauchy-Green  Operator.  The  principle  of 

os  coni 

virtual  work  (2. A)  with  v^  6  C6(ft)  c  E  (ft)  implies  that  the  equilibrium 

fields  u^  are  weak  solutions  of  the  system  of  partial  differential  equa~ 

tions  Cf  (u)  +  F.  ■  0  in  !1,  If  the  body  is  homogeneous,  as  is  assumed 

x  j ,  j  1 

in  this  section,  then  the  system  may  be  written 

(2.14)  »lk  "k  +  F1  -  » 

where 


(2.15) 


Aik  "  cijk£  3  /9xj  3xi 


The  matrix  differential  operator  (Aik)»  with  coefficients  that  satisfy 
the  positivity  and  symmetry  conditions  (1.4),  (1.5),  will  be  called  the 
Cauchy-Green  operator.  Conditions  (1.4),  (1.5)  imply  that  (A^)  is 
strongly  elliptic  (c^^  ^  Cj  J4  0  for  all  non-zero  ni,  and 
hence  elliptic  (det  (c^j^  Sjj)  f1  0  for  all  non-zero  [12,  p.  20], 

G.  Fichera  [5]  has  used  the  theory  of  elliptic  boundary  value  problems 
to  prove  both  interior  and  boundary  regularity  theorems  for  weak  solutions 
of  (2.14).  The  interior  and  boundary  regularity  properties  of  solutions 
wLFE  that  are  implied  by  Fichera' s  results  and  methods  are  described  here 
briefly. 


Interior  Regularity  of  Solutions  wLFE.  Fichera' s  interior 
regularity  theorem  [5,  p.  36]  implies  the  following  results. 

Theorem  2.1.  Let  8  c  RJ  be  an  arbitrary  domain.  Let 
uA  e  Ljnt(n),  eij(u)  €  L*nt(ft)  and  F^  6  L”’int(ft)  where  m  >  0  is  an 
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integer.  Assume  that  (2.4)  holds  for  all  v.^  e  C*(ft).  Then 


u.  S  L^+2,int(a) 


Corollary  2.2.  Let  ft  C  R3  be  an  arbitrary  domain  and  let  u  be  a 
solution  wLFE  of  the  equilibrium  problem  for  ft  with  F^  e  L°,Com(n).  Then 
u±  6  LfH2’int(fl)  . 

Corollary  2.3.  If  the  hypotheses  of  Theorem  2.1  or  Corollary  2.2 
hold  then  u.  €  £?“(«). 

i 

Corollary  2.4.  Let  ft  C  Rs  be  an  arbitrary  domain  and  let 
u  £  E^°c(ft)  satisfy  e.  (u)  -  4  (u.  +  u,  .)  *  0  in  L2°c(ft).  Then  there 

1  j  L  J-  >  J  Jl1 

exist  constants  a^  b^^  such  that  u^x)  *  aj^  +  £ijk  bj  xk  in 

Fichera  proved  Theorem  2.1  in  [5]  under  the  hypotheses  f  6  L®(ft), 
u  €  L2(ft).  However,  the  theorem  as  stated  above  is  n  immediate  conse¬ 
quence  of  his  theorem.  Corollary  2.2  is  a  special  case  of  Theorem  2.1. 
Corollary  2.3  follows  from  Theorem  2.1  and  Sobolev's  imbedding  theorem 
[5,  p.  26].  Corollary  2.4  may  be  verified  by  noting  that  u  is  a  solution 
wLFE  in  ft  with  body  forces  F^  ■  0  in  ft.  Thus  uA  e  C°°(ft)  ,  by  Corollary 

2.3,  and  u  +  u  ■  0  in  ft.  The  proof  that  every  such  u  has  the  form 
i» J  j » i  1 

uA  -  a^^  +  bj  x^  is  classical  [11,  p.  71]. 

Boundary  Regularity  of  Solutions  wLFE.  Fichera' s  theorems  on 
regularity  at  the  boundary  imply  the  following  results  (see  [5,  Chapters 


10  and  12]). 

Theorem  2.5.  Let  ft  C  RJ  be  a  domain  with  boundary  3ft  e  C°*.  Let 
u  be  a  solution  wLFE  of  the  equilibrium  problem  for  ft  with 
F±  6  Lj°m(ft)  n  C "(IT).  Then  ui  e  C°°(ft)  and 

(2.16)  c^u)  nj  “  0  on  3^ 


where  n^  is  the  unit  exterior  normal  field  on  3ft. 


Corollary  2.6.  Let  x0  e  3f2  and  assume  that  there  is  a  neighbor¬ 
hood  N^(x0)  ■  {x  |  |x  -  x0 1  <  6}  such  that  3f2  n  N.(xQ)  €  C°.  Moreover, 
let  F  €  L2°m(fi)  n  c“(JT  n  N6(x0)) .  Then  u±  €  C°°(?2  n  N6(x0))  and 
oij(u)nj  •  0  on  3ft  n  N^Cxj). 

Corollary  2.6  is  an  immediate  consequence  of  Theorem  2.5  since 
boundary  regularity  is  a  local  property.  Boundary  regularity  results 
can  also  be  proved  when  3ft  and  F^  have  a  finite  number  of  derivatives. 
The  following  results  can  be  proved  by  the  methods  of  [5];  see  also  [1] . 

Theorem  2.7.  Let  ft  c  Rs  have  a  boundary  point  xQ  such  that 
k+2 

3ft  n  Ng(xfl)  €  C  for  some  6  >  0  where  It  >  0  is  an  integer.  Let  u  be  a 
solution  wLFE  of  the  equilibrium  problem  for  ft  with  F^  e  L^^ft)  n 
L^Cfl  n  Nfi(x0)).  Then  u±  6  lJ1*©  n  Sfi(x#)). 

Corollary  2.8.  Under  the  hypotheses  of  Theorem  2.7, 
u^  €  Ck(ft  n  N^(xo)).  Moreover,  if  k  >  1  then  ay(u)n^  *  0  on 
3ft  n  N^(x0) . 

Corollary  2.9.  Let  ft  c  R3  be  a  domain  with  boundary  3ft  e  C1^2, 
k  >  0.  Let  u  be  a  solution  wLFE  of  the  equilibrium  problem  for  ft  with 
F^  €  l^» com(fl) .  Then  ui  S  C^(f2) .  Moreover,  if  k  >  2  then  Uj^  is  a 
classical  solution  of  the  equilibrium  boundary  value  problem  with  body 
forces  F^  €  C^  2(ft)  C  L^* COm(Sl)  and  zero  surface  tractions;  i.e.,  uA 
satisfies  (2.10)  and 

C2-17>  cijkl  “k.ji,  +  Fi  ■  0  in 

Bodies  whose  boundary  3ft  is  a  piece-wise  smooth  surface  with 
piece-wise  smooth  edges  with  corners  are  of  great  interest  for  applica¬ 
tions.  A  class  of  bodies  of  this  type  are  the  C-domains,  defined  and 
studied  by  N.  Week  [24].  Solutions  wLFE  in  such  domains  are  regular  and 
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satisfy  the  boundary  condition  (2.16)  near  smooth  points  of  3(2,  by 
Corollary  2.8.  At  edge  and  corner  points  of  3(2  condition  (2.16)  is 
meaningless,  because  n^  is  undefined,  and  the  only  regularity  property 
that  remains  is  the  LFE  condition.  For  this  reason  the  LFE  condition  is 
sometimes  called  the  "edge  condition"  [20]. 

Equilibrium  Problems  with  Kon-Zero  Surface  Tractions.  The  formu¬ 
lation  (2.4)  of  the  principle  of  virtual  work  is  appropriate  for  the  case 
of  zero  surface  tractions.  The  surface  traction  at  a  point  xQ  €  3(2  is 
by  definition  the  vector  o^j(u(x0))  nj(xtt)  and  hence  is  defined  only  at 
boundary  points  where  the  boundary  values  (Jy (u(x0))  and  the  normal  vector 
nj(x0)  exist.  If  a  portion  S  c  3(2  is  sufficiently  smooth  for  n^  and 
boundary  values  of  0^  (u)  to  exist  on  it  then  the  principle  of  virtual 
work  can  be  extended  to  include  the  boundary  condition 


(2.18) 


°ij(u) 


t^  on  S 
0  on  3(2  -  S 


To  do  this  the  term 

(2.19)  -J  t^  v^  dS  *  Work  done  against  surface  tractions 
must  be  added  to  (2.4),  so  that  the  extended  principle  becomes 

(2.20)  J  0^00  (v)  dx  -  J  Fi  v±  dx  -  J  t±  v±  dS  -  0 

for  all  v  €  Ecom((2).  Moreover,  it  is  known  from  Sobolev's  imbedding 
theorem  that  every  v  6  Ecom((2)  has  boundary  values  v  S  L2  (S)  on  smooth 
portions  S  C  3(2  (1,  p.  38],  In  the  important  special  case  where  3(2  Is 
piece-wise  smooth  then  (u)  n^  exists  almost  everywhere  on  3(2  and  S 
may  be  replaced  by  3(2  in  (2.18),  (2.19)  and  (2.20). 
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3.  UNIQUENESS  THEOREMS  FOR  PROBLEMS  WITH  PRESCRIBED  SURFACE  TRACTIONS 

The  strain  energy  theorem  for  classical  solutions  of  the  elasto- 
static  equilibrium  problem  with  body  forces  and  zero  surface  tractions 
states  that  [18,  p.  173] 

(3.D  wn  “  2  jn  aij(u)  eij(u)  dx  "  I  !0  Fi  ui  dx 

The  uniqueness  of  classical  solutions  is  a  corollary.  In  this  section 
the  strain  energy  theorem  is  extended  to  arbitrary  domains  ft  and  all 
solutions  wFE  (*  solutions  wLFE  and  zero  stresses  at  infinity  if  ft  is 
unbounded)  and  solutions  wLFE  and  zero  displacements  at  Infinity.  The 
uniqueness  of  solutions  wLFE  with  prescribed  stresses  or  displacements 
at  infinity  follow  as  corollaries.  The  simple  case  of  solutions  wFE  is 
treated  first. 

Theorem  3.1.  Let  u  be  a  solution  wFE  of  the  equilibrium  problem 
with  body  forces  F^  e  L20m(ft)  and  zero  surface  tractions  in  a  domain 
ft  C  R3.  Then  the  strain  energy  equation  (3.1)  holds. 

The  proof  is  Immediate  from  the  representation  (1.7)  for  and 
the  definition  of  solution  wFE,  since  one  nay  take  ■  ui  e  E(ft)  in 
(2.4). 

Corollary  3.2.  Uniqueness  of  Solutions  wFE.  Let  u^1^,  u£2^  be 
two  solutions  wFE  of  the  equilibrium  problem  with  the  same  body  forces 
F.^  €  Econi(ft)  and  zero  surface  tractions.  Then 

(3.2)  0lj(u<l))  “  °ij(u(2>)  in  n 

and  there  exist  constant  vectors  a^,  b^  such  that 
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(3.3)  u<‘>  (x)  -  u£2^(x)  -  ai  +  e  k  bj  xk  in  ft 

Proof.  -  u<‘>  -  u<2>  is  a  solution  wFE  wirh  body  forces 
F^  =  0  in  ft  and  zero  surface  tractions.  Thus  (3.1)  holds  with  ■  O 
and  o  (u)  ■  0  in  L2(ft)  by  the  positive-definiteness  of  the  energy. 
Moreover,  (u)  G  C°°(ft)  by  Corollary  2.3  and  hence  fl^(u)  (x)  =  0  in  ft 
which  implies  (3.2).  Finally,  Corollary  3.4  implies  u^(x)  ■  a^+G^^b^x^ 
which  implies  (3.3). 

Corollary  3.3.  Uniqueness  of  Solutions  wLFE  with  Prescribed 
Stresses  at  Infinity.  Let  ft  C  RJ  be  unbounded  and  let  u^1^,  u£2^  be  two 
solutions  wLFE  of  the  equilibrium  problem  with  the  same  body  forces  F^, 

CO 

zero  surface  tractions  and  the  same  stresses  at  infinity.  Then  (3.2) 
and  (3.3)  hold. 

Proof .  By  hypothesis,  both  O^(u^)  -  and  (^(u^)  - 
are  in  L2(ftR  ^  for  some  R  >  0.  It  follows  that  the  difference  field 
»  ui^  “  ui^  satisfies  C  (u)  e  L2(ftR  „)  •  Hence  u^  is  a  solution 
wFE  with  body  forces  F^  =  0  in  ft  and  zero  surface  tractions.  Equations 
(3.2),  (3.3)  follow  as  in  the  proof  of  Corollary  3.2. 

The  uniqueness  theorem  for  solutions  wLFE  with  prescribed  dis¬ 
placements  at  infinity  will  be  based  on  the  following  generalization  of 
Theorem  3.1. 

Theorem  3.4.  Let  u  be  a  solution  wLFE  of  the  equilibrium  problem 
with  body  forces  F^  €  Lj0m(ft)  and  zero  surface  tractions  in  an  unbounded 
domain  ft  c  R3 .  Moreover,  let  u  satisfy 

(3.4)  )  !ul"2  dr  -  +• 

Jr  r,° 
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for  some  R  >  0  and  6  >  0.  Then  u  is  a  solution  wFE  in  ft  and  the  strain 
energy  equation  (3.1)  holds. 

A  proof  of  Theorem  3.4  is  given  at  the  end  of  the  section, 
following  the  statement  and  discussion  of  the  remaining  uniqueness 
theorems. 

Corollary  3.5.  Uniqueness  of  Solutions  wLFE  with  Prescribed 
Displacements  at  Infinity.  Let  ft  C  R3  be  unbounded  and  let  u£  ^ ,  uj2^ 
be  two  solutions  wLFE  of  the  equilibrium  problem  with  the  same  body 

QO 

forces  F^,  zero  surface  tractions  and  the  same  displacements  u^  at 

infinity.  Then  (3.2)  and  (3.3)  hold. 

Proof.  By  hypothesis  iu^  -  u°°l_  r  *  0( r1^2) ,  r  “,  k  ■  1,  2. 
—  T  f  0 

It  follows  by  the  triangle  inequality  that  the  difference  field 
u^  «  «<»  -  uf>  satisfies  Iulr  g  -  0{x^2) ,  r  ■+•  “,  or  equivalently 

(3.5)  ,uSr,5  "  0(r)*  r  ■*”  “ 

which  implies  condition  (3.4).  Moreover,  u  is  a  solution  wLFE  with 
Fa  ■  0  and  zero  surface  tractions.  Hence  (3.1)  holds  with  Fj^  «  0,  by 
Theorem  3.4,  and  the  conclusions  (3.2),  (3.3)  follow  as  before. 

Uniqueness  Theorems  for  Problems  with  Non-Zero  Surface  Tractions. 
The  uniqueness  theorems  proved  above  are  valid  for  arbitrary  bounded  and 
unbounded  domains  ft  c  R3 .  No  local  or  global  restrictions  are  imposed 
on  ft  or  3ft.  If  a  portion  S  C  8ft  is  smooth  enough  for  the  surface 
tractions  0^ (u)  nj  and  surface  integrals  (2.19)  to  be  defined  then 
solutions  wLFE  with  non-zero  surface  tractions  t^  on  S  are  defined  by 
the  principle  of  virtual  work.  The  uniqueness  theorems  for  solutions 
with  zero  surface  tractions  extend  immediately  to  this  case  because  the 
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difference  of  two  solutions  with  the  same  surface  tractions  t^  is  a 
solution  with  zero  surface  tractions. 

Other  Growth  Conditions  at  Infinity.  It  is  clear  from  condition 
(3.4)  of  Theorem  3.4  that  condition  (3.5)  is  only  one  sufficient  condition 
for  uniqueness.  Generalizations  are  obtained  by  replacing  (3.5)  by 

(3.6)  |u,r,6  “  r  ® 

where  p(r)  is  a  function  such  that 

(3.7)  I  p(r)"1  dr  -  +» 

jR 

If  (2  is  an  exterior  domain  ({x  |  |x[  >  R}  C  fi  for  R  >  R0)  and  if 
the  body  is  isotropic  as  well  as  homogeneous;  i.e.,  (1.8)  holds,  then 
the  uniqueness  theorem  can  be  proved  under  weaker  growth  restrictions 
than  (3.4).  Indeed,  under  these  conditions  Fichera  [4]  has  proved  that 

(3.8)  ui(x)  *  0(1)  *»  ui(x)  «  Odxp1)  and  *  ^Clxf2) 

M.  E.  Gurtin  and  E.  Sternberg  [10]  have  rederived  this  result  and  proved 
the  complementary  result  that 

(3.9)  o±  (x)  -  0(1)  -  Ui(x)  -  (KM'1)  and  ai;J(x)  -  0(|x|“2) 

Moreover,  these  results  are  based  on  an  expansion  theorem  for  biharmonic 
functions  in  a  neighborhood  of  infinity  and  are  independent  of  30.  Thus 
the  uniqueness  theorems  for  solutions  wLFE  with  prescribed  displacements 
or  stresses  at  infinity  in  homogeneous  isotropic  solids  are  valid  for 
arbitrary  exterior  domains  0  under  the  conditions 


(3.10) 


u^x)  -  uA(x)  -  0(1),  (x| 
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and 

(3.11)  cr^CuHx)  -  o“j(x)  *  0(1),  |x|  -*•  » 
respectively. 

Proof  of  Theorem  3.4.  The  idea  of  the  proof  is  to  put  v^,  • 
in  the  principle  of  virtual  work  identity  (2.4),  as  in  the  proof  of 
Theorem  3.1.  However,  this  cannot  be  done  directly  when  u  is  a  solution 
wLFE  because  v  6  ECOm(fi)  must  have  compact  support.  Instead,  let 
v^(x)  ■  <Kx)  u^x)  where 

(3.12)  <Kx)  -  (K(|x|  -  R)/6) ,  R  >  0,  6  >  0,  x  €  R* 

and  <l>  €  C  (R)  is  a  function  such  that  i|/'(t)  <  0,  0  <  iJ>(t)  <  1  and 

’  1,  T  <  0 

(3.13)  *(T)  -  < 

0,  T  >  1 

These  properties  imply  that  (J>  €  C*(R3)  ,  0  <  <{>(x)  <  1  and 

1,  |x|  <  R 

(3.14)  <p(x)  ■  < 

0,  jxl  >  R  +  6 

It  follows  that  for  all  u  6  E^0C(fi),  v  ■  $u  €  ECom(ft)  and 

(3.15)  »lfJ  -  +  ♦, j  “i 

Moreover, 

(3.16)  ♦ fj(x)  -  #’((|x|  -  R)/S)  x j / 6 | x | 
and 

(3.17)  supp  C  fiR>0(> 


*> 


With  this  choice  of  v. 


(3.18)  e^v)  "  $  eij^u)  +  1  (*ti  uj  +  ui> 


and  hence 


(3.19) 


°ij(u)  eij^  "  ^  °ij^u^  eij^u^  +  °ij^u^  ^,i  uj 

-  <f>  cr^u)  e^j  (u)  +  S'1  i{/'(|x|  -  R)/<5)  a±j(u)  x±  Uj 


where  x^  «  Xj/|x].  By  assumption  F^  e  Lj°m(f)) .  Choose  R0  so  large  that 
supp  C  (x  |  | x |  <  R0J  and  substitute  v^^  *  $  Uj^  and  (3.19)  in  (2.4) 


with  R  >  R„ .  The  result  can  be  written 


J  <J>  0±j(u)  vij(u)  dx  +  6 


(3.20) 


+  6-1  f  r 

J»R,6 

-  jB  ri  "i  * 


o^u)  5^  Uj  dx 


dx  -  0 


The  goal  of  the  remainder  of  the  proof  is  to  calculate  the  limit 
of  equation  (3.20)  for  R  -*■  «  and  to  show  that  the  limiting  form  is  the 
energy  equation  (3.1).  To  this  end  define 

(3.21)  f (R)  -  |  iK5_1(|x|  -  R))  0^ (u)  e±j (u)  dx  -  F±  u±  dx,  R  >  R„ 

By  equation  (3.20)  an  alternative  representation  is 

(3.22)  f  (R)  -  -5-1  |  ifi'(6-l(|x|  -  R))  (u)  x±  dx 

JjJR,5 

The  properties  of  f(R)  that  are  needed  to  complete  the  proof  of  Theorem 
3.4  are  described  by 

Lemma  3.6.  f  e  C^R,,®)  and  has  derivative 

(3.23)  f'(R)  ■  -6~l  |  ij»'  (6”J  ( |x|  -  R))  C^u)  e^u)  dx  >  0 
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In  particular,  f(R)  is  monotone  non-decreasing  on  [R0,«>).  Moreover, 


(3.24)  f  2(R)  <  u0r.6  f'(R).  R  >  R0 

where  M2  -  (6"1  Cj)  1^*  (t>  I  - 

Proof  of  Lemma  3.6.  Form  the  difference  quotient 


h"l{f(R  +  h)  -  f (R) }  -  f  h"1(v(o‘I(|x|-R-h))  -  \pC6~1 ( |x|-R))} 


(3.25) 


R,R+h+6 


x  ^(u)  eij^  dx 


The  quotient 


(3.26)  h‘,{*(«“1(|x|-R-h))-*(«“l(|x|-R))}  -5“x  *' (S’1  ( |x|-R)) ,  h-  0 


uniformly  for  x  in  bounded  sets  in  R3.  Moreover,  (^(u)  ei^(u)  is 
Lebesgue  integrable  on  bounded  subsets  of  fi.  Thus  passage  to  the  limit 
h  ■*  0  in  (3.25)  is  permissible  by  Lebesgue' s  dominated  convergence  theorem. 
Hence  f'(R)  exists  for  all  R  >  R„  and  is  given  by  (3.23).  It  is  easy  to 
show  that  the  integral  in  (3.23)  defines  a  continuous  function  of  R  which 
is  non-negative.  The  monotonicity  of  f(R)  follows. 

To  prove  the  inequality  (3.24)  note  that  (3.22)  implies  the 

estimate 

(3.27)  |f(R)  |  <  6'1  |  |*’(«-I(|x|-R))!  |a1:)(u)  Uj|  dx,  R  >  R„ 

j£2R,6 

Moreover,  by  repeated  application  of  Schwarz's  inequality 

(3.28)  1  CT±j <\i)  Ujj  <  (a±j(u)  xi  <?kj(u)  xk^l/2  <uj  uj^1/2 

f  3  >1/2 

(3.29)  !«„(■)  *jl 
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f 

5  .i1/2r  i 

!  U/2 

°ij(u)  S1  akj(u)  \  z  l 

L  (Vo>  v  J  [J 

■  <Vu)  %>2J 

(3.30) 


3  3 


-  I  (a.,.(u)  til  0^00  "  °...,(u)  a.,(u) 

j-1  J  i-1  j-1  1J  13 

Now  e^  -  together  with  (1.9)  imply 

(3.31)  cx  a±i  0±j  <  e1j  “  ^ijU  °ij  ak£  1  co  °ij  aij 


for  all  0^  ■  °ji'  C°mbining  these  inequalities  gives 
(3.32)  Idi^u)  x±  Uj|  <  cj/2  (^±j(u)  e±j(u))1/2  (u^  u^)^2 


Substituting  in  (3.27)  and  using  Schwarz's  Inequality  again  and  equation 
(3.23)  gives 

If 00  I  <  «“*  cj/2f  |«,(6"1(Ix|-R))|  (oi;j(u)  e1:j(u))1/2  (Uj  Uj)1'2  dx 

%,6 

(3.33)  <  6"1  c1/2  (J  |*' |  Ojj  (u)  e±j(u)  dx]^  fj  |\P' |  UjU^dxj1^ 

^R.d 

<  6_1cj/2yl/2  (6  f  *  CR) ) 2/2  lulR>6 

where  y  «  Max  |ij>'(x)|.  Squaring  (3.33)  gives  (3.24). 

Proof  of  Theorem  3.4  Concluded.  Lemma  3.6  implies  that  f(+®) 
exists  as  a  finite  number  or  +°°.  It  will  be  shown  that  f(-H“)  *  0.  There 
are  three  cases  to  consider. 

Case  1.  0  <  f(+°°)  <  +00.  In  this  case  there  exists  R.  >  R.  such 
‘  ~  •  ~  —  1  —  0  • 

that  f(R)  >  f(Ri)  >  0  for  R  >  Rt .  Hence  (3.24)  can  be  written 

«•»>  'St  (ffe)  ■Fwi""  •  *i*i 

and  integration  gives 


v 
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«-35>  fib'  -  ffe  -  I,  ,u;Mdr  ■  *  1  *' 

R1 

In  particular,  since  f(R)  >  0  for  R  >  Rj , 

M2  fR  -2 

(3.36)  >  j  lu|  zs  dr  for  R  >  Rx 

But  this  contradicts  hypothesis  (3.4)  of  the  theorem.  Hence  Case  1 
cannot  occur. 

Case  2.  f(+°°)  <  0  and  f(Rj)  “  0  for  some  Rx  >  R0 .  In  this  case 

0  <  f (Rj)  <  f(4®)  <  0;  i.e.  f(-H»)  -  0. 

Case  3.  f(+“)  <  0  and  f(R)  <  0  for  all  R  >  R„ .  In  this  case 
(3.34)  and  (3.35)  hold  and  the  latter  can  be  written,  since  |f(R)|  *  -f(R), 

<3-373  jmr 5  iwr +  m"!  C  '<s dr- *  i  *> 

Hence  condition  (3.4)  implies  that  1(4®)  ■  0. 

It  has  been  shown  that  (3.4)  implies  f(+“0  *  0;  that  Is, 

(3.38)  lira  f  t|»(6_1(|x|-R))  a  (u)  e  (u)  dx  -  [  F.  u.  dx 

R-H~  h  13  13  h  1  1 

Since  ip(6-J ( [x|  -  R))  is  a  monotone  increasing  function  of  R  for  each 
fixed  x  €  R3  and  tends  to  1  everywhere  when  R  ®,  (3.38)  implies  equation 
(3.1).  In  particular  <  °°  because  F^  u^  dx  is  finite.  This 
completes  the  proof . 


4.  UNIQUENESS  THEOREMS  FOR  OTHER  EQUILIBRIUM  PROBLEMS 


The  purpose  of  this  section  is  to  show  how  the  methods  and  results 
developed  above  can  be  extended  to  the  most  general  equilibrium  problems 
of  linear  elastostatics.  Equilibria  subject  to  other  boundary  conditions, 
equilibria  in  inhomogeneous  anisotropic  bodies  and  n-dimensional  general¬ 
izations  are  discussed.  In  each  case  the  boundary  conditions  for  dis¬ 
placement  fields  wFE  and  wLFE  are  defined  by  appropriate  subspaces  of 
E(ft)  and  E  (ft) ,  respectively,  and  a  corresponding  form  of  the  principle 
of  virtual  work  is  given.  Regularity  and  uniqueness  results  for  the  new 
problems  are  indicated  without  proofs.  In  fact,  the  proofs  of  sections 
2  and  3  are  valid  with  minor  modifications. 

Equilibrium  Problems  with  Prescribed  Surface  Displacements.  The 
case  of  zero  surface  displacements  is  discussed  first.  Suitable  subspaces 
of  displacements  fields  are 

(4.1)  Efl(ft)  ■  Closure  in  E(ft)  of  Ecom(ft)  n  {u  |  supp  u  c  ft) 

(4.2)  E g°C(ft)  “  Closure  in  E^OC(ft)  of  ECom(ft)  n  {u  |  supp  u  C  ft) 

The  topologies  in  E(ft)  and  E^0C(ft)  are  those  defined  by  (1.28)  and  (1.29), 
respectively.  The  notation 

(4.3)  E0Com(ft)  -  ECOm(ft)  n  E#(ft) 

is  also  used.  A  solution  wFE  of  the  equilibrium  problem  with  body  forces 
6  Lj°m(ft)  and  zero  surface  displacements  in  a  field  u  e  E0(fl)  that 
satisfies  (2.4)  for  all  v  6  E0 (ft) .  Similarly,  a  solution  wLFE  of  the 
same  problem  is  a  field  u  S  Eg0C(ft)  that  satisfies  (2.4)  for  all 


35 


v  e  Eg0m(ft).  Problems  with  non-zero  surface  displacements 


(a. 4)  u^x)  -  fi(x)  ,  x  e  an 

may  be  reduced  to  the  preceding  problem  if  there  exists  a  field 

u°  €  E^°C(ft)  n  {u  |  (u°)  £  Lz0tn(ft)}.  Then  u^  -  ui  -  *s  a  solution 

wLFE  with  zero  boundary  displacements. 

The  remaining  boundary  conditions  can  be  formulated  only  when  3ft 
is  piecewise  smooth.  It  will  be  assumed  that  8ft  is  a  C-domain  in  the 
sense  of  [24].  For  such  domains  the  unit  exterior  normal  field  n^(x)  is 
defined  and  continuous  at  all  points  of  3ft  except  edges  and  corners  and 
one  can  define  the  normal  and  tangential  components  of  vector  field  on 


3ft  by 

(4.5) 

ui  "  ui  +  ui  »  ui  “  (uj  vi 

Moreover, 

V  T 

*  0  for  all  u^,  v^  and  hence 

(4.6) 

V  V  .  T  T 

Ui  vi  “  ui  vi  +  ui  vi 

Equilibrium  Problems  with  Prescribed  Tangential  Surface  Tractions 
and  Normal  Surface  Displacements.  Suitable  subspaces  of  displacement 
fields  are  defined  by 

(4.7)  Ev(ft)  -  E(ft)  n  {u  |  uV  -  0  on  3ft} 

(4.8)  E J°c(ft)  -  E^°c(ft)  n  {u  |  uV  -  0  on  3ft} 

The  existence  of  uV  and  uT  on  3ft  for  all  u  €  E^oc(ft)  follows  from  Korn's 

inequality  and  Sobolev's  imbedding  theorem.  A  solution  wFE  of  the 

coin  y- . 

equilibrium  problem  with  body  forces  F^  €  Lz  (ft),  zero  tangential 
surface  tractions  and  zero  normal  surface  displacements  is  a  field 


*> 
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u€  Ey((2)  that  satisfies  (2.4)  for  all  v  e  £^((2).  Similarly,  a  solution 
wLFE  of  the  same  problem  is  a  field  u£  (H)  such  that  (2.4)  holds  for 
all  v  €  E^°m((2)  •  Ev((2)  n  ECOn(^).  Problems  with  non-zero  surface  trac¬ 
tions  and  displacements  are  treated  by  reducing  them  to  the  preceding 
case  through  subtraction  of  a  suitable  field. 

Equilibrium  Problems  with  Prescribed  Normal  Surface  Tractions 
and  Tangential  Surface  Displacements.  This  problem  is  dual  to  the 
preceding  one.  Appropriate  classes  of  displacements  are 

(4.9)  ET(fi)  -  E((2)  n  {u  |  uT  ■  0  on  8(2} 

(4.10)  E*oc((2)  -  Ei'oc(n)  n  {u  |  uT  *  0  on  8(2} 

Equilibrium  Problems  with  Elastically  Supported  Surface.  Physi¬ 
cally,  this  corresponds  to  the  case  where  surface  displacements  produce 
surface  tractions  that  satisfy  Hooke's  law: 

(4.11)  a^(u)  n^  +  3  ui  «  0  on 

JLoc 

where  3  >  0  is  defined  on  8(2.  A  solution  wLFE  is  a  field  u  6  E  ((2) 
such  that 

(4.12)  [  0. . (u)  e,.(v)  dx  -  I  F,  v.  dx  +  6  u.  v.  dS  ■  0 

'(2  ij  h  1  '8(2  1  1 

for  all  vS  Ecom(S2) .  Identity  (4.12)  is  the  principle  of  virtual  work 
for  this  problem,  the  last  term  being  the  virtual  work  done  against  the 
induced  surface  tractions  by  the  virtual  displacement  v.  It  follows  from 

(4.12)  that  (4.11)  holds  at  smooth  points  of  8(2. 

Equilibrium  Problems  with  Nixed  Boundary  Conditions.  A  mixed 
problem  that  includes  the  preceding  problems  as  special  cases  can  be 
formulated  by  decomposing  8(2  into  five  portions  and  imposing  one  of  the 
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boundary  conditions  defined  above  on  each  portion.  Thus,  if 

(4.13)  3JI  »  Sj  U  S2  u  Sj  u  u  S5  (disjoint  union) 
and 

(4.14)  E^°C(Q)  =  E£oc(H)  n  (u  |  u  =  0  on  Sj ,  uV  -  0  on  S2 ,  uT  =  0  on  S3} 
then  the  principle  of  virtual  work 

(4.15)  I  0. .  (u)  e. .(v)  dx  -  [  F.  v.  dx  +  B  u.  v.  dS  *  0 

Jn  ij  Jn  1  •>s5  1  X 

fo.  all  v  €  Ecom(£))  ■  E£oC(Q)  n  ECom(fi)  characterizes  the  solutions  of 
m  m 

the  equilibrium  problem  that  satisfy  u  ■  0  on  S1 ,  uV  “  0  and 

(o^  (u)  Hj)T  «  0  on  S2,  uT  =  0  and  (o^  (u)  n.)V  =  0  on  S3,  (u)  n^  ■  0 

on  S4  and  0^  (u)  +  8  =  0  on  S5 . 

Regularity  and  uniqueness  theorems  will  be  discussed  for  this 
mixed  problem  since  it  includes  the  others  as  special  cases. 

Regularity  Theorems.  The  interior  regularity  properties  of 
solutions  wLFE  of  the  mixed  problem  follow  from  Theorem  2.1  and  are 
exactly  the  same  as  for  the  case  discussed  in  section  2.  Concerning 
boundary  regularity,  it  can  be  shown  by  the  methods  of  Fichera's 
monograph  [5]  that  if  ft  is  a  C-domain  of  class  C°°  such  that  S®  -  interior 
of  in  is  a  C°°  manifold  for  k  -  1,***,5,  and  if  F.^  e  C°°(ft)  n  l2°m(n) 
then  solutions  wLFZ  of  the  mixed  problem  satisfy 
Ui  6  c"(n  u  s[  u  S’  u  s|  u  sj  u  S®)  n  L2’£oc(fl).  The  condition 
u^  6  L2,£°C(ft),  which  follows  from  Korn's  inequality  and  Sobolev’s  therre 
is  the  "edge  condition"  that  is  needed  for  uniqueness.  The  boundarv 
conditions  on  S2 ,  S3  and  S$  are  not  discussed  by  Fichera  in  [5]  b.: 
be  treated  by  his  methods. 
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Uniqueness  Theorems.  Solutions  vFE  of  the  mixed  problem  lie  In 

(4.16)  En(fl)  ■  E(ft)  n  {u  |  u  -  0  on  ,  uV  ■  0  on  SJt  uT  ■  0  on.  S,} 
end  satisfy  (4.15)  for  all  v  €  E  (ft).  The  strain  energy  theorem  for  the 

WL 

problem  is 

(4.17)  \  Jq  0±iM  e±j(u)  dx  +  \  ^  dS  -  f  F±  u±  dx 

vhere  the  first  equation  defines  the  strain  energy  for  the  mixed  problem. 
The  uniqueness  of  solutions  vFE  is  an  Immediate  corollary.  Solutions 
with  prescribed  stresses  or  displacements  at  Infinity  will  be  defined  by 
(2.9)  and  (2.10),  respectively,  as  In  the  surface  tractions  problem. 
Moreover,  the  strain  energy  theorem.  Theorem  3.4,  extends  to  solutions 
wLFE  of  the  mixed  problem.  In  fact,  the  same  proof  Is  valid  because  If 
u  €  E^°C(ft)  and  4  €  c"(R*)  then  v  -  $u  6  E^°“(ft)  -  Em(ft)  n  Econ(ft) .  The 
uniqueness  of  solutions  vLFE  of  the  mixed  problem  with  prescribed  dis¬ 
placements  at  infinity  Is  an  Immediate  corollary.  It  can  also  be  shown 
that  the  displacement  fields  for  the  mixed  problem  are  unique  except  In 
the  special  case  of  the  pure  surface  tractions  boundary  condition  (S„ « 3ft) , 
Inhomogeneous  Bodies.  The  uniqueness  and  energy  theorems  given 
above  remain  valid  If  the  constant  stress  strain  tensor  Is  replaced 

by  a  field  c^^Cx)  that  Is  Lebesgue  measurable  in  ft  and  satisfies  (1.9). 
The  Interior  and  boundary  regularity  theorems  of  section  2  are  valid 
when  Cjju^x)  has  sufficient  differentiability  in  ft  and  ft,  respectively; 
cf.  [I,  p.  132J. 

n-Dlmenslonal  Problems.  Flchera  [5]  has  developed  his  theory  for 
an  n-dlmenslonal  generalization  of  the  equations  of  elastostatlcs.  All  of 
the  theorems  given  above  extend  to  this  n-dlmenslonal  problem  with  only 


rotational  changes.  The  cases  n  ■  1  and  n  -  2  are  applicable  to  elasto- 
static  fields  that  are  functions  of  only  one  or  two  of  the  Cartesian 
coordinates . 


5.  A  DISCUSSION  OF  RELATED  LITERATURE 


Fichera's  paper  [4]  of  1950  provided  the  first  significant  exten¬ 
sion  of  Rirchhoff 's  uniqueness  theorem  to  unbounded  domains.  His  result 
(3.8)  implies  that  equilibrium  fields  in  homogeneous  isotropic  bodies  in 
exterior  domains  have  finite  energy  if  the  displacements  vanish  at 
infinity.  The  uniqueness  of  equilibrium  fields  in  such  bodies  is  an 
lmnediate  corollary.  Corresponding  results  for  fields  vhose  stresses 
vanish  at  infinity  follow  from  the  1961  result  (3.9)  of  Curtin  and 
Sternberg  [10] .  The  author  knows  of  no  genaral  uniqueness  results  for 
anisotropic  bodies  in  exterior  domains  or  for  bodies  vhose  boundary  is 
unbounded. 

In  Fichera's  monograph  [5]  of  1965  the  existence  and  uniqueness 
of  classical  solutlona  to  elastostatic  equilibrium  problems  in  bounded 
domains  with  smooth  boundaries  is  proved  by  the  methods  of  functional 
analysis.  This  provides  an  alternative  to  the  classical  Integral  equa¬ 
tion  methods  cited  in  the  introduction.  However,  the  formulation  and 
techniques  employed  by  Flchera  can  provide  more  general  results. 

Fichera's  semi-weak  solutions  (Lecture  7)  are  essentially  the  solutions 
wFE  of  this  paper.  Hence,  Fichera's  results  (Lectures  7  and  12)  Imply 
the  uniqueness  of  solutions  wFE  for  bounded  domains  and  boundary  condi¬ 
tions  for  which  Korn's  inequality  is  valid.  For  the  zero  surface  dis¬ 
placements  problem  the  inequality  holds  for  every  bounded  domain.  For 
the  zero  surface  tractions  problem  it  bolds  for  domains  with  the  cone 
property. 
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The  literature  on  uniqueness  theorems  in  linear  elastostatics  up 
to  1970  was  surveyed  in  a  monograph  by  R.  J.  Knops  and  L.  E.  Payne  [13] 
published  in  1971.  This  work  also  contains  uniqueness  theorems  for  a 
class  of  weak  solutions.  However,  the  hypothesis  that  the  displacement 
fields  are  continuous  in  H  restricts  the  scope  of  these  results. 

Uniqueness  theorems  for  plane  crack  problems  were  proved  by 
J.  K.  Knowles  and  T.  A.  Pucik  in  1973  [14]  under  the  assumption  that  the 
displacements  are  bounded,  but  not  necessarily  continuous,  at  the  crack 
tips.  The  elegant  differential  Inequality  method  used  in  this  work 
provided  the  inspiration  for  the  proof  of  Theorem  3.4. 

The  methods  employed  in  this  paper  to  prove  uniqueness  theorems 
for  solutions  wLFE  in  arbitrary  domains  were  introduced  by  the  author 
during  the  period  1962-64  in  a  series  of  papers  on  boundary  value 
problems  of  the  theory  of  wave  propagation  [25,  26,  27].  The  article 
[27]  contains  as  a  special  case  uniqueness  theorems  for  elastodynamic 
problems  in  arbitrary  domains. 
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COMPLETENESS  OF  THE  EIGEN  FUNCTIONS  FOR 
GRIFFITH  CRACKS  IN  PLATES  OF  FINITE  THICKNESS 


Introduction. 

E.S.  Folias  [1]  has  constructed  the  displacement  and  stress  fields 
near  a  Griffith  crack  as  an  expansion  in  eigenfunctions.  Hie  eigen¬ 
functions  were  derived  by  an  operational  method  due  to  Lur€  [2]  and  the 
question  of  completeness  arises.  The  purpose  of  this  report  is  to 
prove  the  completeness  by  a  constructive  method.  The  method  employed 
is  to  solve  the  boundary  value  problem  by  Fourier  analysis  and  to 
evaluate  the  resulting  integrals  as  residue  series.  The  terns  in  these 
series  are  precisely  the  eigenfunctions  used  by  Folias. 


J 
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Notation. 

A  system  of  Cartesian  coordinates  (x,y,z)  is  used, 
occupies  the  region  defined  by 

-°°  <  X  <  ®  ,  -®  <  y  <  OO  t  -h<Z<h 

The  crack  is  defined  by 

-c  <  x  <  c  ,  y  3  0  ,  -h  <  z  <  h 


_ _ _ 


The  plate 


2 


The  components  of  the  displacement  field  in  the  stressed  plate  are 

u(x,y,z)  ,  v(x,y,z)  ,  w(x,y,z) 

The  corresponding  stress  tensor  components  are 


_  _  .  r3u  3v  9Wv  ,r  3u 
°x  "  X(57  +  3y  +  3Z3  +  2G  3x 


_  „  rr3u  .  8v.  T 

Txy  G^3y  3x^  Tyx 


T  -  G(|H  +  f£)  -  T 
XZ  '•dz  3xJ  zx 


0  -x(|h  +  |y:  +  |W)  +  2G|v 
y  '■Sx  3y  3zJ  3y 


T  =  *  t 

Tyz  Ul3z  3yJ  zy 


_  .  ^  f3u  3v  3w>  7r  3w 
°z  X(57  +  37  3I3  2G  31 


The  displacement  field  satisfies  the  field  equations 


32u 

K  32,J  H 

k  32u 
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3 

,3u  . 

3v 

J? 

k  — 7  H 

ay 

I  + 
3z 

a 

57 

C57  + 

3y 

3ZV 

32v 

„  32v 

2 

3 

,3u  . 

3v 

~ 7  H 
3x^ 

— 7  H 
3y 

.i? 

a 

3y 

(57  + 

3y 

32w 

.  32w  , 

t  32w 
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,311  . 

3v 

J? 

k  — J  * 

3y 

“7  + 
3z 

a 

dZ 

te  * 

ay 

-  0 
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where  a  -  .  The  boundary  conditions  are 
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Txz  “  0  »  ryz  *  0  »  °z  “  0  for  <  x,y  <  “  »  z  -  ±h 

V  0  •  V ' 0  ’  °y  “  "°0  for  -c  <  x  < c  •  y  - 0i  •  1*1  4  •> 

u,v,w  -  0(1)  for  x2  *  y2  +  «  ,  |z|  <  h 

Symmetries . 

u(x,y,z) 
v(x,y,z) 
w(x,y,z) 

If  follows  from  the  symmetries  that 

v ' 0  •  It  " 0  ■  It"  0  f°r  1*1 >  c  >  y  - °  • 

for  |x|  >  c  ,  y  =  0  , 

for  |x|  >  c  ,  y  *  0  ,  |z|  <  h 


o  iSL-  0 

3y  u  »  3y  u 


whence 


r  «  0  .  t  =0 
xy  *  zy 


|z|  <h 
|z|  <  h 


X 

y 

z 

odd 

even 

even 

even 

odd 

even 

even 

even 

odd 

Moreover,  if 
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then 


u1 (x,0,z)  =  ^  u(x,+e,z) 

[u] (x,0,z)  =  u+(x,0,z)  -  u' 

(x,0,z) 

=  [|^]  =  0  for  -«»  <  x  <  »  ,  | 

z|  <  h 

*  =  0  ^or  -00  <  x  <  “  ,  | 

z|  <  h 

for  -»  <  x  <  00  ,  | 

z|  <  h 

etc. 


even  in  y 


while 

[0]  =  2(0)+  ,  [v]  =  2v+  ,  etc.}  odd  in  y 

Note:  these  vanish  at  points  of  continuity. 

Application  of  the  Fourier  Transform  in  x  and  y  . 

Define 


u(p,y,z)  =  - — 1/2 

(2  *)1/Z  i 


e'ipx  u(x,y,z)dx 


etc. 


w 

0(p,q,z)  -  -±in  }  e-iqy 


z)dy 


^  |  |  uCx.y.^dxdy 


etc. 


Then 


However, 
the  crack. 

f  c  c2(ig 


and  hence 
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(* h' 


ipu 


3X 


2- 
-p  u 


'•ay-'  3y 


etc. 


i,v,w  and  their  derivatives  may  have  discontinuities  across 
Note  that,  if  £^(y)  e  L^QQ  for  k  *  0,1,2  and 
n  C2(R_)  and  f(0±)  ,  f 1  (0±)  are  finite 


00  0 

(f’Mq)  -  -\/2  (  f  ♦  [ 

(2*) V  n  1 


)  e-1^  f  •  (y)dy 
0- 


i-s/2  (e"iqy  f(y)]  +  e'iqy  f(y)]_  +  iq  f  e‘i(^  f(y)dy) 

!Tr)i/z  0+  * 


(2tt) 


=  -  y/2  +  ^ 

(2t01/2 


(f")"(q)  -  -  ^\/2  +  iq  (f'f  (q) 

'  ‘  w1'2 ' iq  ^ ' q2  kq) 


These  results  and  the  symmetries  (p.  3)  imply 


I 


6 


,3u.“  .  ~ 

(gy)  =  iq  U 


,»  *  ,2.1/2,3u.A  2  ~  ,,3u.A 

'  V  ‘^0  '  q  U  t(^o 


<|jff  (P.0*,«» 
etc. 


,3v.~  „  ,2.1/2  A  .  ~ 

<5y>  ©  v0  +  iq  v 


,32v.~  .  ,2.1/2  A  2  ~ 

^T>  *  -  iq  (f )  vQ  -  q  v 


0"  -  iq  w 


,32w.~  _  ,2.1/2,3w  A  2  ~ 

tl*  ‘  <5*  (3y>  ‘  q  w 


3y 


Taking  the  Fourier  transform  of  the  field  equations  (p.  2)  and  using  the 
above  results  gives 


>2. r. 


u  2  ~  2  ~  ,2.1/2,3u>A  2,  2  -  ~  ,2.1/2  A  .  dw. 

d  7  -  P  u  ‘  *  u  -  ^  +  a  (-  P  u  -  pq  v  -  ip(f)  vQ  +  xp  jj.)  -  0 


,2~ 

d  v 


dz 


-  p2  v  -  q2  v  -  iq(|)1/2  vQ  +  a2(-  pq  u  -  q2  v  -  iq(|)1/2  vQ  +  iq  ^)  *  0 


d2w  2  ~  2  ~  ,2.1/2 ,3w.A  ,  2..  du  .  .  dv  ,2. 1/2 ,3V..1  .  d2w.  _  n 

■£z  -  p  w  -  q  w  -  <?>  +  a  dp  a? +  iq  a?  “  W  ^  0 


or 


1 

i 


7 


2~ 


d 

dz 


u  ,  2  2»~  2  ,  ~  ~  .  diV\  ,2«2/2r  ✓Su-i*  .  2  ~  , 

j  -  (p  +  q  )u  -  a  p(pu  +  qv  -  i  jj)  -  C^)  '  [(^)  +  ip  a  vQ] 


,2~ 


TT'  (P2  +  "  a2q(pu  +  qv  -  i  ^)  -  (|)1/2[iq  vQ  iq  a2  vfl] 


dz 


,2- 


a*a2)  ^ -  (pV> *i4|*q|)-  cf)1/2[cg)"  * a2<H>j 


3yy 


Note  that 


f(x)  is  even  «-*■  £(p)  is  even  and 


£(P)  =  (|)1/2  |  cos  px  f(x)dz  -  Fcf(p) 
0 


w 

f(x)  =  (|)1/Z  |  cos  px  fCp)dp 


while 


f(x)  is  odd  f(p)  is  odd  and 


f (P)  “  -i  (|)V2  j  sin  px  f(x)dz  -  -i  Fsf(p) 


f(x)  -  i  C |)V2  j  sin  px  f(p)dp  -  (|)1/2  |  sin  px  Fsf(p)dp 


The  analogous  formulas  are  valid  for  functions  of  y  .  Hence  the 
symmetries,  p.  3,  imply 


-i 


8 


u  *  -i  U  ,  v  *  -i  V  ,  w  ■  W 
where  U  ,  V  ,  W  are  real-valued.  In  fact. 


IM  W 

U  *•  ”  I  j  sin  px  cos  qy  u(x,y,z)dxdy 
0  0 

CO  00 

V  ■  ^  f  |  cos  px  sin  qy  v(x,y,z)dxdy 
0  0 


W  ■  —  I  I  cos  px  cos  qy  w(x,y,z)dxdy 
0  0 

Hence  the  differential  equations  for  u  ,  v  ,  w  on  p.  7  are  equivalent 
to 


d^ 

37 


a2p  ^  -  a2p(p  U  +  qV)  -  (p2  +  q2)U  -  c|d1/2[FsC^o  -  pa^] 


-0  -  a2q  g  -  a2q(p  U  ♦  q  V)  -  (p2  ♦  q2)V  -  (|)1/2[  -  (l+a^qF^] 

n<tfl2,  dU  *  n  dV.  -2  .  n2.w  .  ,2.1/2rF  ,3w.  2,  ,3v.  , 

(l+a  )  j  +  a  (P  +  q  )  '  (P  +  q)W-(f)  [^(57)  +  »  Fc^n ■ • 


dz 


This  can  be  written  as  a  2nd  order  3x3  matrix  system  of  ODE's, 
namely 


[ul 

U 

fFll 

L 

V 

-  (p2  ♦  q2) 

V 

■ 

F2 

w 

w 

W 

where,  if  H  ■  (U,V,W)r  (T  *  transpose) 


ur-  +  bS*  cct 

d7  3* 


fl 

0 


0  0 

1  0 

0  l+a' 


B  - 


[a  p  a  q 


0  -a2p' 

0  -a2q 
2 


-aV 

-a2pq 


-a2pq  0 
-a2q2  0 


Note  tliat 


AT  -  A  ,  BT  -  -B  ,  CT  -  C 

It  follows  that  L  is  formally  self ad joint  with  respect  to  the  scalar 
product 

*2 

(0,V)  ■  J  0^  V  dz 

n 


In  fact,  integration  by  parts  gives 


(7.L7) 


where 


If  the 


J  (O7  A  7'  +  CT  B  7'  ♦  UT  C  7)dz 


Z2  *2 


UT  A  7'  ♦  UT  B  V)  -  J  (U»T  a  7*  ♦  U,T  B  7  - 


-1  2 


•  UT  A  V'  -  U  'T  A  V  ♦  UT  B  7] 


L 

J  (tf'T  a  -  rr,T  b  ♦  uT  c)  7  dz 


[U,7]  ♦  f  (aD"  +  b  U1  ♦  C  D)t  7  dz 

z-,  l 


-1  z 


[U,7]  +  (L  U,  7) 

z-, 


[7,7]  -  UT  A  7’  -  U,T  A  7  +  UT  B  7 

index  notation 

U  -  OIl,U2,U5)T  .  7  -  (V^V^Vj)7 


lT  C  7)dz 


is  used  the  bilinear  form  [7,7]  can  be  written 
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[U,V]  -  UXV'  +  u2V£  +  Cl+a2)U3V^ 

-  uivi  -  U2V2  -  C1*«2)UJV3 

-  a2pUxV3  -  a2qU2V3  +  a2pU3Vx  ♦  a2qU3V2 

Boundary  Conditions  Associated  with  L  . 

The  symmetry  properties  of  the  displacement  field  wrt  z  (p.  3) 
imply  that 

tf^z-O  "  0  *  ^z-O  "  0  *  Wz»0  *  0 

It  follows  that 

*  o  ,  -  0  ,  V(0)  “  0 

Note  that 

B.C.l  1^(0)  -  0  ,  U£(0)  -  0  ,  U3(0)  -  0 

is  selfadjoint  for  L  ;  i.e. 

U  and  7  satisfy  B.C.l  •*  •  0 

The  B.C.’s  at  z  *  ±h  imply  corresponding  B.C. ’s  for  U  .  To  write 
then  note  that  (p.  6) 


since 


UU 


V  ■  G4  *  i<?  5)  ■  -i  G(g  -  I  <0 


~  , ..  .  ,2.. 1/2  *  dw..  _  _  dw 

c  -  X(ip  u  +  lq  v  •  (?)  v0  +  3z5  2G  Hz 


X(p  U  +  q  V)  +  (X  +  2G)  -  (|)1/2  X  v( 


»Ip  U  *  q  v  .  «jli  gj  -  (£)1/2  X  ^ 

a  -I 


X-  “  ,  *-.-2  , 

m-  z  A  A  ^ 


2  »  ,  ma2  -  2a2  ■  m  ,  m(a2  -  1)  *  2a2 


_  »  2a 

n  ~2 —  »  25 
a  -1 


X  a2-l  m  _  2a2 


*  m- 


a  +i 


It  follows  that 


-  p  W(h)  -  0  ,  -  q  W(h)  -  0 


P  uo.)  ♦  <,  voo  ♦  ^  -  t|)i/2  ;0cM 

a  -1 


Note  that 


B.C.2  UJ(h)  -  p  U3(h)  -  0  ,  U£(h)  -  q  U3(h)  *  0 


pU,(h)  +  q  U,(h)  +  ^  "  0 

1  z  a  -1  5 


is  also  selfadjoint  for  L  .  In  fact,  if  U  and  V  satisfy  B.C.2  then 


[U,V]2=h  -  p  U:V3  +  q  U2V3  -  (aZ-l)U3  Cp  \  *  q  V2) 


-  p  UZV1  -  q  U3V2  +  CaZ-l)  Cp  ^  +  q  U2)V3 


a2p  U2V3  -  a2q  U2V3  +  a2p  +  a2q  U3V2  =  0 


BV  Problem  for  U  -  (U^ ,U2 ,1 


,W)T  . 


where 


L  U  -  (p2  +  q2)U  »  F  ,  0  <  z  <  h 


MqUCO)  ♦  N0U’  CO)  -  0 
MjU(h)  +  NjU' Ch)  -  CTCh) 


'  FsOo  ’  pa2  Fc  v0 
c|)1/Z  -  Cl+a2)  q  Fc  v0 

-  Fc0o  -  a2FcC|l)0 


0  0  0 


10  0 


Mq  -  0  0  o  ,  N0  -  0  1  0 

0  0  1  l  0  0  0 
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0  0  -p 

O 

O 

«h- 

0  0  -q 

’  V 

0  10 

7 

p  q  0 

0  0  a*+l 

* 

k 

a2-l 

and 


CM  - 


$1/2fc  v0 


Method  of  Solving  the  BV  Problem  for  U,V,W  . 

—  2  2  —  — 

To  solve  the  BV  problem  the  general  solution  of  L  U  -  (p  +  q)U**F 
and  BC  at  z  *  0  will  be  constructed  as  a  function  of  the  parameters 


“o  -  UCO)  ,  v0  -  V(0)  ,  wj  =  ^21 

The  B.C.'s  at  z  ■  h  will  then  be  used  to  calculate  u^  ,  ,  w^  . 

Solutions  of  the  Equations  L  U  -  (p  +  q  )U  =  0  . 

This  equation,  written  in  terms  of  components  (U^t^U^)  ■  (U,V,W)  , 
is  obtained  from  the  system  on  p.8  by  setting  the  right-hand  side  equal 
to  0  .  Note  that  this  system  coincides  with  Lurg,  p.  150  (3.2.123 
under  the  correspondence 


-i  U  —  u  ,  -i  V  «■*  v  ,  \i~  w 

ip  —  ,  iq  —  82  ,  -  (p2  +  q2)  —  D2 

ip(-i  U)  ♦  iq(-i  V)^'PUMV^-*0 


Usr4  has  given  a  complete  solution  of  the  system  L  U  -  (p  +  q  )U  in 
equation  (3.2.15),  (3.2.17).  To  adapt  them  to  the  present  notation  write 

S2  -  P2  +  q2  ,  s  -/ pZ  +  q2  >  0  ,  D  *  is 

cos  2D  ■  cos  isz  =  cosh  .;z  ,  sin  zD  *  sin  isz  =  i  sinh  sz 

Then  the  solution  (3.2.15)  becomes 

2 

U  =  (cosh  sz)u0  +  2  -  (p2u0  +  pqvQ  +  pw*) 


V  -  (cosh  sz)vfl  +  \  2  -S^-—  (pqu0  ♦  q2v0  +  qw’) 


2  .. 


— ■  wq  +  X  ~S~  '  z  cosh  sz)  (Pun  +  3vn  +  WA) 


0  M'0  O'* 


This  solution  satisfies  B.C.l  at  z  ■  0  .  The  solution  (3.2.17)  becomes 


u  -  3i^sz  +  _a  (Lsirftsz  +  ^coshsz)(p2  +  +  ps2 


2 (a  +1)  -i  s” 


v  -  pLy'  ♦  -V-  )Cpqu.  +  q2v-  ♦  qs2wn) 


2 (a  +1)  -s‘ 


W  »  (cosh  sz)w0  -  (pu,  +  qv.  ♦  s2w0) 


This  solution  satisfies 


B.C.l' 


U(0)  -  0  ,  V(0)  =  0  ,  -  0 


Solution  Basis  for  L  U  «  (p^  +  q^)U  . 
uQ  -  1  ,  v0  =  w0  =  0  gives 


Similarly 

Uq  ■  1  ,  v^-0,  Wq  ■  0  gives 


(z  cosh  sz  -  IHiLlijL 
s  2(a  +1)  s  sZ 


_4_  (z  CoSh  sz  - 
2(a  +1)  s 

a2  z  sinh  sz  _ 

- 2 - c - P 

2(a%l)  s 


un  =  0  >  v'  =  i  ,  wn  =  0  gives 


’  u5 ' 

Vs 

8 

w5 

—S, -  (z  cosh  sz  -  lML2i)Ea 

2(a  .1)  s  s 

iJSLII  *  -4—  (z  cosh  sz  - 

s  2 (a  +1)  s  sZ 


a  z  sinh  sz 


u0  *  v0  “  0  »  wq  =  1  8ives 


a(a  +1) 


a  (a  +1) 


,  ,  sinh  sz-. 

Cz  cosh  sz - - - )  p 


a_  (z  cosh  sz  -  SHiLSi)  q 
2(a  *1)  s 


I  Yr  I  I  cosh  sz - n - sz  sinh  sz  i 

2(az+l) 

It  is  evident  from  the  B.C.  at  z  =  0  that  these  six  solutions  are 
linearly  independent  and  hence  span  the  solution  space  of 
LH«  (p2  +  q2)U  . 

Solutions  of  L  U  -  (p2  +  Q2)U  =  F(z)  . 


Solutions  of  L  U  -  (p  +  q  jU  =  F(z)  . 

The  variation  of  constants  formula  will  be  used.  For  this  purpose 


1 


it  is  convenient  to  write  the  equation  as  a  1st  order  system.  The 
equation  has  the  form 


A  U"  +  B  tJ’  -  CsU  =  F 


where 


C  =  s*l  -  C  = 
s 


r  2  2  2  2 

1  s  +  a  p  a  pq 


2  2  2  2 
a  pq  s  +  a  q 


Now  A  =  A  >  1  ,  whence 


A  =  A1/2  A1/2  ,  A1/2  =  (A1/2)T  = 


f  1  0  0  1 

0  10 

0  0 


Thus 


A1/2  U"  +  A-1/2  B  U'  -  A'1/2  C  U  =  A'1/2  F 


Put 


r 


V  = 


U, 


U 


2 


U  =  A 


Then 


V"  + 


(A‘1/2  B  A'1/2)Vf 


or 


where 


V*'  +  ba  V'  -  cA  V 


Ba  *  A‘1/2  B  A'1/2  = 
CA  *  A_1'/2  C3  A  1/2  = 


Z  *  A 


■1/2  F 


Explicitly, 


ba- 


(l+a2)1/2 


0 
0 
l  P 


0 

0 

q 


2  ,  2  2 
s  +  a  p 


a2pq 


2 

a  pq 


2  2  2 
s  +  a  q 


0 


0 
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A  1st  order  system  equivalent  to  the  above  2nd  order  system  may  be 
obtained  by  setting 

Y  *  7  ,  Z  =  V' 

X  =  $  -  (7,) 

Then 


Y1  »  F'  «  Z  ,  z*  «  V"  -  -  B^Z  ♦  C.Y  +  E 

and 


0 

1  ' 

o' 

X'  = 

X  + 

•CA 

“ba 

s 

or 


X’  -  M  X  +  H(z)  ,  H  »  (^) 

where 


M  » 


A  fundamental  matrix  for  X'  =  H  X  is  a  6  x  6  matrix  solution  4>(z) 
of 


’  4-*  (z)  -  M  $(z) 

.  «K0)  -  1 


An  explicit  representation  of  0(z)  can  be  derived  from  the  solution 
basis  for  L  U  ■  s%  .  Indeed,  each  solution  tP  ■  (iP  of 


- 
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L  tP  ”  s2lP  gives  a  solution  of  X^ '  *  M  X-*  ,  namely 


■yj  ' 

fA1/2  Jfi 

Vj* 

s 

A1/2  IP  * 

Thus,  in  view  of  the  B.C.'s  at  z  *  0  , 


$(z)  = 

CX1  X2  (1+a2) 

'1/2X6  X4  Xs  Cl+aV^X3) 

M2 

”l 

U1 

u2 

(l+a2)‘1/2U6 

U4 

u5 

(l+aV^V 

V1 

V2 

C  r1/2v6 

V4 

V5 

(  )'1/2V3 

(l+aV^V 

c  )VV 

w6 

(l+a2)1/2W4 

(l+a2)1/2^ 

W3 

U1' 

u2. 

( r1/2u6* 

u4* 

U5' 

c  r^V* 

V1* 

V2’ 

(  3-1/2/ » 

v4' 

/» 

(  )‘1/V* 

(l+a2)^</2Wlf 

c 

W6* 

( )VV* 

(  )1/2  IV5* 

W3' 

where  iP  ,  V3  ,  Vp  are  defined  on  pp.  16-17. 

The  fundamental  matrix  makes  it  possible  to  calculate  a  solution  of 
X'  =  M  X  +  H  ,  namely 

z 

X(z)  -  $(z)’J  CO  H(C)dc 
0 

Indeed, 

z 

X’  -  i'1^)  Hfcjdc  +  «(z)  *_1(z)  H(z) 

0 

■MX  +  H 


- 


Moreover, 
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X(0)  -  0 


Thus  the  general  solution  of  X*  *  M  X  +  H  is  given  by 

z 

X(z)  -  $(z)  XQ  +  *(z)  J  KC)'1  H(c)d£ 

0 

The  direct  calculation  of  is  difficult,  but  note  that  if 


10  -1J 


and 


then  X'  -  M  X  - 

E’(z) 


E(z)  -  j  XCz)T  P  X(z) 


X(z)T  P  X'(z)  *  X(z)T  PM  XCz)  =  0 


because 


and  hence 


(PM)T 


-  PM 


Thus  E(z)  =*  const.  ¥•  solutions  of  X’  *  M  X  .  Take  X(z)  ■  *(z)  Xq 
(Xq  e  R®  arbitrary) .  Then 

2E(z)  -  (*(z)  X0)T  P  *(z)  Xq  -  Xj  *(z)T  P  *(z)  Xq 

«  xj  P  XQ  »  2E(0)  V-  X0  e  R6 
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It  follows  that 

*(z)T  P  *(z)  -  P  ¥  z  e  R 
Since  P  is  non-singular 

(P"1  *(z)T  P)  $(z)  »  1 

whence 

♦(z)"1  =  P_1  *(z)TP 

Thus 

z 

(*)  X(z)  -  *(z)  Xq  ♦  *(z)P-1  J  *(C)T  P  H(Odc 

0 

Solution  of  the  B.V,  Problem  of  p.  13. 

Recall  that 

X  ■  $,)  •  (Ux  U2  (l*a2)1/2  U3  UJ  UJ  Cl*a2)1/2  UJ)T 

*0  ’  <u0  v0  »**2)1/2  «0  4  vj  (l*a2)1/2  wj) 

Thus  the  solution  (*)  satisfies  B.C.l  (at  z  *  0)  u^  *  -  wQ  »  0  . 

Thus  if  we  write 

z 

X*(z)  -  *(z)  P'1  |  $(C)T  P  H(t)dt 
0 

then  substituting  in  (*)  gives 
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Ux(z)  -  *n(z)  uQ  ♦  *12(z)  Vq  ♦  (l+a2)1/2  *16COw’  ♦  xj(z) 

U2(z)  -  *n(z)  uQ  ♦  422(z)  vq  ♦  (l+a2)1/2  *26(i)  wj  ♦  Xj(z) 

(l+a2)1/2U3(z)  -  431(z)  uq  +  *32(z)  Vq  ♦  (l*a2)1/2  *36(z)  w*  ♦  xj(z) 

or  (see  p.  21)  if  Uu(z)  -  xj(z)  ,  V*(z)  -  x£(z)  ,  Wff(z)  -  (H-a2)"1/2X*(z)  , 
U(z)  -  U1 (z)  uQ  ♦  U2(z)  v0  +  U3 Cz)  ♦  UK(z) 

V(z)  -  V*(z)  Uq  ♦  V2(Z)  Vq  +  V*(z)  W*  +  V*  (z) 

W(2)  -  ^(z)  Uq  +  W*(z)  Vq  ♦  (z)  W*  ♦  W* (z) 

Thus  (p.  15) 

2 

U(z)  *  (cosh  sz)  Uq  +  Z  -Siy  —  (P?U0  +  WVq  +  pw^)  +  lf(z) 

V(z)  -  (cosh  SZ)  Vq  ♦  ^  (pqu0  ♦  q2v0  +  qw')  -  V*(z) 

W(z)  -  S-^-SZ  «o  +  V  ^  SZ  '  z  cosh  sz)  (pu0  +  qv0  +  WP  *  ^ (z) 

To  complete  the  solution  of  the  B.V.  problem  of  p.  13  the  initial  values 
u0,v0*w0  must  c^osen  so  t^at  B.C.2  at  z  *  h  is  satisfied.  The 

derivatives  Uf,V*,W*  are  needed.  They  are  given  by 

U*  (z)  -  s(sinh  sz)Uq  +  £  (sinh-SZ  +  sz,  cosh  si }  ^  +  ^  +  ^  +  ^.(z) 

V’ (z)  -  s(sinh  sz)vQ  +  ^  (sinh  52  *-SZ  C°*h  SZ)(pqu0  ♦  q2vQ  ♦  qw*)  +  Vw*(z) 

W*  (z)  -  (cosh  sz)Wq  +  %£■  (  -sz  sinh  sz)(puQ  +  qvQ  +  w£)  ♦  ^’(z) 

Thus  the  B.C.2  (p.  13)  gives 
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U’  (h)  -  p  W(h) 

2 

,  .  ,  ,  a  ,sinh  sh  +  sh  cosh  sh%  ,_2 .  _  A  _.ts 

*  s(sinh  sh)ufl  +  -j-  ( - g - - )  (p  uQ  +  pqvQ  +  pw'Q) 

-  p  sin^-sh  v’Q  -  (S-1T1|----  -  h  cosh  sh) (puQ  +  qvQ  +  wfl) 

+  if'  (h)  -  p  Vf  (h)  =  0 

V'  (h)  -  q  W(h) 

2 

.  s(SM  sh)v0  .  (pqUg  *  q2v„  *  qwj) 

.  q  -  S|i  -  h  cosh  shJtpu,,  *  qv0  ♦  wj) 

+  V^'GO  -  q  l<f(h)  -  0 

p  U(h)  +  q  V(h)  +  W  (h) 
a -1 

-  p(cosh  sh)uQ  +  E|_  h  sing-s^  (p2u0  +  pqvQ  +  pw’) 

+  q(cosh  sh)vQ  +  ^2“  h~SilV"~  ^u0  +  qZv0  +  qw0^ 

2  2,2 

+  (^±)  (cosh  sh)w^  +  C^1)  -f  C-sh  sinh  sh)  (puQ  +  qvQ  +  w') 

^  2  A. 

+  p  lf(h)  +  q  V*(h)  *  C^)  ^'OO  «  (|)1/2  v0(h) 

a  -1 

This  is  a  system  of  linear  equations  for  uq>vq»wq  of  the  form 
dllu0  +  d12vQ  +  c^jWq  -  f^(p»q)  (M  -U1Tr  (h)  +  p  vf  (h) ) 
d2iU0  +  d22v0  +  d23w^  -  f2(p,q)(-  -V^'fh)  +  q  W* (h)) 

2 

d31U0  +  d32V0  +  d33W0  "  f3^P*q)C“  “P  U*  00  -  q  Vff(h)  -  C^W*'  (h) 

*  (|)1/2^0(h)) 


•  # 
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where 


2  2  2  2 

dll“  s(sinh  sh)  +  -  (sinh  sh  +  sh  cosh  sh)  — (sinh  sh  -  sh  cosh  sh) 

2  2 

*  s(sinh  sh)  +  a  p  h  cosh  sh 


dj2  pq  ^sinh  sh  +  ^  cosh  pq  (Sinlj_sh  -  h  cosh  sh) 


a  pq  h  cosh  sh 


l13m  T  p  ^  S"  +  h  cosh  sh)  -  |  ^  sh  -  ^  p  (?-^-  sh  -  h  cosh  sh) 
*  a2ph  cosh  sh  -  p  — 


,sinh  sh 


21  -  V  W  ♦  h  cosh  sh)  -  pq  C 


sinh  sh 


-  h  cosh  sh) 


a  pqh  cosh  sh  -  d12 


a  _2  ,sinh  sh 


a  2  ,sinh  sh 


d22  =  s(sinh  sh)  +  C  ™  +  h  cosh  sh)  "  V  q  C~s  -  h  cosh  sh) 


2  2 

*  s(sinh  sh)  +  a  q  h  cosh  sh 


2  2 
j  _  a  _  ,sinh  sh,  , _ t  „u->  _  sinh  sh  a  ,sinh  sh  t _ i_ 

d23  ~T  q  ( — s -  “  cosh  sh'  ”  q  — I - T  q  £ — S - h  cosh  SW 


-  sinh  sh  +  a2qh  cosh  sh 


d31  -  p  cosh  sh  +  p3  j  sinh  sh  +  ^-  pq2  j  sinh  sh  -  sh  p  sinh  sh 

a  -1 


p  cosh  sh  +  (i-  psh  -  (— «t— -)  i -  psh)  sinh  sh 
L  a -1  1 

2  2  2 
p  cosh  sh  +  -U  psh  (  -  —n — )  sinh  sh  *  p  cosh  sh  — 7 —  p  sh  sinh  sh 
L  a  -1  a  -1 


,a2+l,  a2 


'  9 


2  2 

j  _  a  _2_  h  sinh  sh  _ _ u  „u  .  a  .3  h  sin 

dJ2  «  fj-  p  q - g - ♦  q  cosh  sh  +  -j-  q  - g 

2  2 
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2  2  2 
i  _3  h  sinh  sh  ,a  +1,  a  , 

q  - - - (_ — )  sh  q  sixth  sh 

L  5  a-1  L 


a2  a2  a2+l 

■  q  cosh  sh  +  4-  qsh  sinh  sh  -  (■■■---)  qsh  sinh  sh 

1  1  a-1 

a2 

-  q  cosh  sh  +  -■*-  qsh  sinh  sh  ) 

L  a-1 

a2 

*  q  cosh  sh  — n —  qsh  sinh  sh 
a-1 

.  £  p2  hji^jh  ,  £  ,2  hsighsh  .  a^l  C05h  sh  .  ^1,  ^  ta  sinh  sh 

2  s  L  s  a-1  a-1  2 

7 


2  2  2  -  2 

-  V  sh  sinh  sh  -  %-  sh  sinh  sh  +  cosh  sh 

2  2  a-1  a2-l 

“2  a2+l 

+  cosh  sh 


2  2 
=  (-^-)  sh  sinh  sh  + 
a-1  L  a-1 

2  2 
-  i-li  cosh  sh  — % —  sh  sinh  sh 
a-1  a-1 


The  Cofactors  of  Q  »  (d^J  . 

Let  Qjk  *  (cofQ)jjj 

Q11  "  d22d33  "  d32d23 

2  2 

-  (s  sinh  sh  +  a2q\  cosh  sh)  (- A  cosh  sh  -  -4 —  sh  sinh  s) 

a-1  a-1 

2 

♦  (q  cosh  sh  -  -4 —  q  sh  sinh  s)  (+  sinh  sh  -  a  qh  cosh  sh) 

a-1  s 

»  ■<}—  [(a2+l)s  sinh  cosh  -  a2s2h  sinh2  +  (a2+l)a2q2h  cosh2  -  a*q2h2s  sinh  cos 
a-1 

2  4 

♦  q2[~  sinh  cosh  -  a?h  cosh2  — % —  h  sinh2  +-§—  h2s  sinh  cosh] 

s  a-1  a-1 

\  .... 


-  C(a>HaX--  -  aVh)  cosh2  ♦  ( -XpL  -  sinh: 

a -X  a-1  a-1 


*  s  -  .  4  .  ii&i)  sinh  cosh 

a-1  a-1  s  a-1 


"1 
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.  2  cosh2  -  xffia2?!?.  sioh2  *  -  PZ  sihh  cosh 

a  -1  a -1  s(a  -1) 

-  ♦  ?A2h-  a2Is^2;h  cosh2  +  j^sVjjLE.2  sirfl  cosh 

a  -1  a -1  s(a -1) 

2,  2  2,,  2  2,  -  2,  2  2s  2 

-  Ul*9Jk  -  2-0.  cosh2  sh  +  2  S-fe  *3  P  sinh  2sh 

a -1  a -1  s(a  -1) 

Q12  =  "  (d21d33  "  d31d235  *  d31d 


l31d235  *  d31d23  "  d21d33 


(p  cosh  -  psh  sinh)  (-  &  sinh  +  a2qh  cosh) 
a -1  s 

2 

—  cosh  - 

a‘-l 


-  (a2pqh  cosh)  (rj-^-  cosh  — | —  sh  sinh) 
a  -1 


a  -x  a  -l 

^  sinh  cosh  +  a2pqh  cosh2  +  5-92*1  sinh2  -  5-P9§*L_  gi^h  cosh 
s  (aM)  a2-l 

-  a— pqh  Cosh2  +  5.  glifr  sinh  cosh 

a  -1  a  -1 

-  .  «!p!>  .  (a2pqh  .  J^S!>  .  cosh2 

a  -1  a  -1  a  -1 

+  (  -  ^)  sinh  cosh 


2  , 

Lpl  4 

a2pqh 

(1  ♦  -y- 

-  iflL)  cosh' 

a -1 

a  -1 

a  -1 

a2pqh 

cosh2sh  - 

^3.  sinh  2sh 

a  -1 

a2-l 

21d32  ' 

d31d22 

2  a 

»  (a  pqh  cosh)  (q  cosh  — j —  q  sh  sinh) 

a  -1 

a2  2  2. 

-  (p  cosh  — n —  psh  sinh)  (s  sinh  +  a  qTi  cosh) 

a  -1 
a2 

*  (cosh  -  -* —  sh  sinh)  (-ps  sinh) 
a -1 

2 


_  *5  0 

a  phsz  siitfi  sh  -  ps  sinh  sh  cosh  sh 


a  -1 
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2  2 

-4 —  phs2  +  -4 —  phs2  cosh2  sh  -  2^  sinh  2sh 
a-1  a-1 


Q21  =  '  (d12d33  "  d32d135  =  d32d13  "  d12d33 
2 

=  (q  cosh  — % —  qsh  sinh)  (a2ph  cosh  -  p 

a-1  s 

2  2 

-  (a2pqh  cosh)  cosh  -  -4 —  sh  sinh) 
a-1  a  -1 


a2pqh  cosh2  -  23.  sinh  cosh  -  ^  ■  sinh  cosh  +  —  sinh2 

s  a-1  a4-l 


.  cosh2  +  sinh  cosh 

a  -1  a  -1 

.2  ...  „  .2_  ,  .2,_2. 


a 

a 


J3!>  .  (a2pqh  ♦  SJ3!l  -  cosh2 

r-1  a-1  a-1 


+  (-  23)  sinh  cosh 
s 

-  -  2ja!i .  ijs!i  cosh2  sh 

a  -1  a  -1 


sinh  2sh  *  Q. 
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3  dlld33  "  d31d13 
=  (s  sinh 

-  (p  cosh  -  a 


Q22  *  alla33  "  u3l' 

inh  +  cosh)  (44^  cosh  -  -4 —  sh  sinh) 

a-1  a  -1 

2  2  m 

cosh  — n —  psh  sinh)  (a  ph  cosh  -  2.  sinh) 
a-1  s 

£il  sinh  cosh  -  4^  cosh2 

'-1  az-l  a-1 


a 


*  s  sinh  cosh  - 

a-1  a-1  -  - 

2  4  2  2  2  2. 

-  a2p\  cosh2  +  2L  sinh  cosh  +  -  2.-  -  ■  sinh  cosh  -  a-& - 
s  a-1  a-1 

.  .  l2.(sV)hCcosh2  .  a2p2h)  cosh2 

a-1  a-1 

2  2 

♦  }X  -  E_  )  sinh  cosh 
a-1  s 

.  .W)ji,  ,  a2„  .  p2  .  sV,  cosh2 

aZ-l  az-l  az-l 


sinh  cosh 


i 


sinh  cosh 


.  sW*»  <■  pV-1) 

s(a-l) 

.  a2CpW)h  .  aVhcosh2  sh  *  jVa&ia!  Sinh 
a-l  a -1  2s  (a-l) 


'  (dlld32  ’  d31d12)  *  d31d12  “  dlld32 


a 


2  2. 


(cosh  -  —  sh  sinh)  (a  p  qh  cosh) 

a' -1 

q  (s  sinh  +  a^p^h  cosh)  (cosh  — | —  sh  sinh) 

a-l 

a2 

(cosh  -  —  sh  sinh)  (-  qs  sinh) 


a-l 


2.2. 


-  qs 


sinh  cosh  +  a  *  sinh2 


a  -1 


.  a2s2qh  +  ffife  cosh2  sh  -  sinh  2sh 
a-l  a-l 


“  d12d23  '  d22d13 

■  (aZpqh  cosh) (-  ^  sinh  +  a  qh  cosh) 

-  (s  sinh  +  a2q2h  cosh)(a2ph  cosh  -  ^  sinh) 

2  2 

-  a  phs  sinh  cosh  +  p  sinh 

2  v 

■  -p  +  p  cosh  sh  -  —  sinh  2sh 


“  "  (d12d23  ‘  d21d13)  “  d21d13  '^d12d23 
■  (a2pqh  cosh)  (a2ph  cosh  -  ^  sinh) 

-  (aZpqh  cosh) (  |  sinh  ♦  a2qh  cosh) 
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»  (a2pqh  cosh) (a^h(p-q)  cosh  +  sinh) 

*  (a2pqh  cosh  sh) (p-q) (a\  cosh  -  j  sinh) 

Q33  =  *11*22  "  *21*12 

-  (s  sinh  +  a2p^h  cosh) (s  sinh  +  a2q2h  cosh)  -  a4p2q2h2  cosh2 

“  s2  sinh2  +  (a2p2sh  +  a2q2sh)  sinh  cosh 
2  2  2  3. 

*  s  sinh  sh  +  a  s"n  sinh  sh  cosh  sh 


|Q|  *  detfd^)  can  be  calculated  from  Lure,  p.  1S3  and  the  correspondence 
(see  p.  14) 


This  gives 


ip  «  »  iq  a2  »  is  ^  D 

iQl  *  2a2h(is)3  sin(ish) (1  +  ^Jj^) 
«  2a 2h  s3  sinh  sh  (1  +  ^^j~) 


Solution  of  the  System  on  p.  25.  We  can  solve  by  means  of  the  relations. 


Thus 
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(*) 


f  IQI  u0  =  ^llfl  +  ^21£2  +  ^31f3 

1^1  v0=  ^12fl  +  ^22f2  +  Q32f3 

l  iQl  w0  13  Vl  +  ^23f2  +  ^33f3 


The  only  real  zero  of  |Q(s)  |  is  at  s  =  /p^+q2  =  0  .  Thus 
V  real  (p,q)  }  (0,0) 

3  Qi5(p,q)f4(p»q) 

uo'p*«  -  £  — 

3  Q,*2(p»q)fiCp,q) 

- 

3  Qi,(p,q)fi(p,q) 

^  — 

Substituting  for  u0,vQ,w'  in  the  equations  on  p.  24  gives 

U(p,q,z)  ,  V(p,q,z)  ,  hr(p,q,z) 

A  A  A 

Residue  Series  Representation  for  u(p,y,z)  ,  v(p,y,z)  ,  w(p,y,z)  . 
The  equations  on  pp.  4-8  give 

00  CD 

u(p,y,z)  =  -1  |  eiyq  U(p,q,z)dq  -  — — 1/2  f  eiyq  U(p,q,z)dq 

(2t01/Z  .1  (2ir)1/Z 

00 

V(p,y,z)  -  -  "1T7?  f  eiyq  v(P»q*z)dq 

C2w r'1  1 


w(p»y,z)  -  -  1  rrj  [  elyq  W(p,q,z)dq 
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In  particular. 


u(p,q,0) 


v(p,y,0) 


^177  |  eW  “o'P'^ 


^17?  1  elyq  votP>")4> 


/V 

3w(p,y,0) 

3z 


wq (p>q)^q 


Now  the  equations  on  p.  32  and  p.  24  give  U,V,W,uQ,v0,w^  as  meromorphic 

functions  of  q  for  each  fixed  p  .  Thus  residue  series  for  the  above 

functions  can  be  obtained  by  deforming  the  contour  in  the  upper  half  of 

the  q-plane  for  y  >  0  (lower  half  for  y  <  0) .  The  poles  of  the 

Integrals  U, . . .  ,w'Q  are  the  zeros  of  |Q(s)  j  .  The  cofactors  Qjk(p,q) 

are  holomorphic  in  the  q-plane.  Examinations  of  the  formulas  for 

u\v\w*  and  f.(p,q)  shows  that  these  functions  are  analytic  everywhere 

-1  -2 

except  at  s  =■  0  ,  because  P  =  0(s  )  .  Thus  special  care  is 
necessary  in  calculating  the  residue  at  s  =  0  (q  **  ±i|p|)  . 


Zeros  of  |Q(s) 1  . 

There  are  two  families  of  zeros 


1)  sinh  sh  =*  -i  sin(ish)  a  0  ish  =  isnh  =  ror  ,  n  *  0,1,2  , 


Thus 


H.  I  ror  .  nu  2  ^  2  fm*2 
s„  ■  *  *q„  •  JK  ■  -1  T  •  P  *  %  ■  -  (-F’ 

qj;  -  -  Cpz  *  cfo2)  ■  % ' 1  •  i 


q0  ■  i  IpI 
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These  are  single  zeros  of  |Q(s)|  for  n  >  1  .  However 

|Q(s)|  »  OCs4)  ,  s-M) 

(P2  +  q2)2  -  (q  *  i|p|)2Cq  +  i|p|)2  ~  (2i|p|)2(q 


Thus  qQ  is,  in  general,  a  higher-order  pole. 

^  ,  .  sin  Zsh  _  ,  ..  sin  2ish  _  .  ,2ish  +  sin  2ish^ 

2)  1  +  - JSK -  =l-i  2ih - 1  C - 2iH - 0 

—  2ish  =  2isvh  =  2Bvh  =  -i3v 

p2  *  <;  -  -  sj  ,  qj  -  -  (p2  ♦  e2)  .  %-  i 


Calculation  of  F(z)  . 

F(z)  is  defined  on  p.  13.  Now  on  y  =  0± 


Thus 


rlHi  =  ^(x,0*,z)  m  _  r3v\ 
0 


3y 


FsC^0 


Fs^0  "  *  P  Fc  T0 


v  riJli  ■  -  F  f— ) 
rcl3y'o  rc'-3zJ0 


-  ilPl)2 


Thus 


-  p  (1-a2)  Fcv0 


(£)1/:1F2Cz)  »  -q  (1+a2)  Fcvq 
(£)1/2F3(z)  =  -  (1-a2)  Fc  (JJ)0 


Calculation  of  if.vW.lf'  .V^'.W77’ 
From  PP.  18-19. 


Write 


Hence 


H(z) 


■  0  ' 

■  0 

[m\ 

a‘1/2Fc*). 

PH(t) 


o 

_ 

’  0  ' 

•  0  • 

[0  -lj 

,-CCz). 

Then 


Hz) 


*Ttt) 


■U  (2) 

»12(z)' 

-21(*) 

»22t2), 

UT« 

♦2lT(« 

12Ta) 

*22T« 

♦'(O  p  n(c) 


*11T  *21Tl 

«12T  ,22T 


'  • 

0 

m  - 

♦21TB 

,-s 

*22TS 

£_1  *TPH 


•c™i 

la 

0 

f*21Tfl 

a  - 

rc:1*21Tm 

A 

o 

-1 

*22T<j 

-*22TG 

*(z)  P'1  *TC0  P  H(0 

f$n(z)  *12(z)j  f-  C‘V1T(0S(0' 

"  U21(z)  *22(z)J  *22Tftm 


,-$11(z)C"1$21TCO^(0  ♦  $12(z)*22T(OSCO' 

+  «22(z)$22T(c)^(0. 


Thus  (pp.  23-24) 

tf( z)  -  (U1T(z),Vir(z),Wir(z))T  =  A‘1/2(XJ,^,X5T 

-  J  A'1/2(*l2(z)*22T(c)  -  *11Cz)C"1t2lT(0)A-1/2F(0dC 
0 

A-VV'cz)  -  (UW)A-1/2 
A'1/2$12(z)  =  (1j4u5D^)a'1/2 

*21TCO  -  CA1/2(u1,u2*tJ6,)A';l/2)T  =  a‘1/2(U1,u2,;D6,)ta1/2 

*22T(C)  -  (A1/2(U4’U5,U3,)A1/2)T  *  A"1/2(D4'U5,U3,)TA1/2 
A~1/2912(z)*22T(t)A~1/2  =  CU4(z)U5(z)U':3(z))A'1CU4,  COU5'  (t)u3t  (C))T 


a-1/2#U(z)C-1#21t(0a-V2  .  (tJ1(z)D2(z)D6(z))Cj1(D1'  (C)U2'  (C)tJ6,  (0)T 
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M^z)  *  (04(z)U5(z)U3(z))  ,  M^z)  =  (01(z)iJ2(z)U6(z)) 

Then 

z 

ifcz)  *  |  {^(zjA^M'CO7  -  ^(zDC^CO^FCOdc 
0 

tf'(z)  =  {M^Cz)A_1M^Cz)T  -  M2(z)C‘1M^(zDT}  F(z) 
z 

*  |  {M*(z)A‘1M’(0T  -  M»(Z)C^(C)T}  F(c)dc 

Similarly, 

tf'(z)  -  A'1/2(xJx*Xg)T 

-  j  A‘1/2(922(z)«22T(0  -  «21(z)C;1021T(c))A"1/2  F(Od; 
A‘1/V2(z)*22TU)A’1/2  =  (04.  (z)U5*  (z)U3’  (z^A^CU4*  COU5'  (C)D3,  (0)T 

A'1/24*21COC^1*21T(i;)A'1/2  -  (U1'  MU2'  COU6*  (Z))A‘1/2C^V1/2CU1'  (OO2,  (;)D6'  (c))' 
z 

if’Cz)  -  j  {MJ(z)A'1M'T(0  -  mjwc^m^co)  FCOdc 

o 

if"(z)  *  {Mj(z)A-1M|(z)T  -  M»(z)Cg1M^(z)T}  F(z) 
z 

+  [  <i^(z)A"IM[Tco  -  Mg(z)c^ft)T)  F(c)dc 
o 

A  If"  (z)  ♦  Btf’Cz)  -  Cslf(z)  -  A{M’(z)A~1M'(z)T  -  M^CjHf^z)7}  F(z) 

-  1 
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An  alternative  derivation  is  as  follows. 

Try  If  (2)  *  ^(2)^(2)  +  ^(2)^(2) 

if’(z)  »  M^(z)r2(z)  +  M^(z)Cj(z) 

+  M1(z)q(z)  +  M^zjqcz)  <-  set  -  0 
If"  (2)  -  M^(Z)U2(2)  ♦  M^(z)q(z) 

♦  M*(z)q(Z)  +  M£(2)q(z) 


A  Tf"  +  B  If '  -  Cs  If  -  A  M[(z)<?*  (z)  ♦  A  M*(z)q(z)  -  F(z) 
Thus 


M^zjqcz)  +  M^Z^'U)  -  I) 
M*(Z)C'(2)  +  M*(z)C*(2)  =  A-1  F(z) 


’u1  u2  u6  u4  u5  u3 ' 

C1 

’  O' 

u2!  u6'  D4'  U5'  D3* 

m. 

A'1? 

Note  that 


*°(z) 


L1 

•1 

l  2J 

1 

< 

♦  (z)  ■ 


A1/2 

0  ' 

n 

-A-V2 

0  1 

.  o 

AV2 

*°(z) 

0 

> 

1 

I-* 
N 
.  »* 

-  C^M£T(z)  F(z) 

A_1M^T(z)  F(z) 
z 

if  (z)  -  J  {Mj^ C2)A_1M*TC0  -  m2(z)C;1m*t(0}  F(Odc 
o 


r^(z)  M*T(z)l  fcs  Oj  fM2(z)  M^z) 
M^(z)  M£T(z) J  (o  -Aj  [m»  (z)  M*(z) 

A  to 

mJ  M|TJ  l-AM^  -AM* 

^CsM2  "  ”2^2  ^  '  M^TAM*' 

-  M{TAM£  mJc^  -  M{TAM* 


Calculation  of  Coefficients  in  the  Residue  Series. 
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1 


[u1  uz  u1 


fu 4  u5  U31 


^(P.q.z)  -  V1  V2  V6  ,  Mjfo.q.z)  -  V4  V5  V3 
W1  W3  W6  IV4  w5  w3 


The  Poles 


%,-id2^*.  sn'-lan’  n  "  1»2»3»*”C°n  “ 

cosh  s  h  =  cos  is  h  *  cosh  ah»  cos  mr  «  f-l)n 
n  n  n  v  J 


sinh  s  h  ■  -i  sin  is  h  ■  -i  sin  mr  «  0 
n  n 


M2(P>qn,h) 


(-l)n 

2(az+l) 

(-»“ 

2  (az+l) 


.  i^l(-i)np  -  5^-1)“^ 

[  aV  ,_„n 


(-l)n 


2(aSi)(-ap 


2  (a  +1)  (-a*) 


M^P.q^h) 


a2hpqn  n  aV(-l)n 

— 2 - — 2 - T~  ® 

2(aSl)(-a‘)  2(az+l)(-a‘J 


^L(-l)n 


M^CP.^.h)  - 


2  (az+l) 
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MiCp,Vh^ 


^P(-Dn  ] 
^C-D" 
(-l)n 


For  a  simple  pole  at  q  “  qg  “  |pl 

Res  {elyq  un(p,q)>  -  Jj®  {Cq-qQ)  e1/q  uQ(p,q)} 
q0  ^0 

-  e"y^  lin  {(q-qQ)  uQCp,q)} 
q^lg 

For  a  double  pole 

Res  (e1™  u0(p,q)>  -  £  (CQ-V2  e1^  u0(p,q)) 

qo 

-  feiyq  ^  «q-V2  uo}  +  iy  ®iyq  {(q_qo)2  “o}] 

-  eiy^  Res  uQ (p,q)  ♦  iy  e'y'p'  a_2(u0) 

q0 

For  higher  order  poles,  correspondingly  higher  order  powers  of  y  appear. 
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